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Abstract 

We study the zero-temperature partition function of the Potts antiferromagnet (i.e., 
the chromatic polynomial) on a torus using a transfer-matrix approach. We consider 
square- and triangular-lattice strips with fixed width L, arbitrary length N, and fully 
periodic boundary conditions. On the mathematical side, we obtain exact expressions 
for the chromatic polynomial of widths L = 5, 6, 7 for the square and triangular lattices. 
On the physical side, we obtain the exact phase diagrams for these strips of width L 
and infinite length, and from these results we extract useful information about the 
infinite-volume phase diagram of this model: in particular, the number and position of 
the different phases. 

Key Words: Chromatic polynomial; antiferromagnetic Potts model; triangular lattice 
square lattice; transfer matrix; Fortuin-Kasteleyn representation; Beraha numbers; confor 
mal field theory. 



1 Introduction 



The two-dimensional (2D) g-state Potts model [1] is one of the most studied models in 
Statistical Mechanics. Despite many efforts over more than 50 years, its exact free energy 
and phase diagram are still unknown. The ferromagnetic regime of the Potts model is the 
best understood case: exact (albeit not always rigorous) results have been obtained for the 
ferromagnetic-paramagnetic phase transition temperature T c (q) (at least for several regular 
lattices), the order of the transition (continuous for < q < 4), the phase diagram, and the 
characterization in terms of conformal field theory (CFT) of the universality classes. (See 
e.g., Ref. [2].) 

The antiferromagnetic regime is less understood, partly because universality is not ex- 
pected to hold in general (in contrast with the ferromagnetic regime). We know the exact 
free energy along some curves of the phase diagram (q, T) (where T is the temperature) for 
some regular lattices [2]; but in this regime, this (partial) solubility of the model does not 
imply criticality (as it is been observed for the ferromagnetic regime with < q < 4). 

The standard g-state Potts model at temperature T can be defined on any undirected 
finite graph G = (V, E) as follows: on every vertex of the graph x G V, we place a spin 
a x taking an integer value in the set {1, 2, ... , q}, where q G N. The spins interact via the 
Hamiltonian 

^Potts(M) = -j J2 ( L1 ) 

where the sum is over all edges e = (x, y) G E joining nearest-neighbor sites, J is the coupling 
constant, and 5 is the Kronecker delta. The ferromagnetic (resp. antiferromagnetic) regime 
corresponds to J > (resp. J < 0). The partition function of the Potts model on a graph 
G at inverse temperature (3 = 1/T > can then be written as 

Z G (q;J(3) = J2^ mPMs (1-2) 
M 

A very useful representation of the Potts model is due to Fortuin and Kasteleyn [3] . They 
showed that the partition function (1.2) can be rewritten as 

Z G (T,v) = « HE ' } vlE ' 1 (1-3) 

E'CE 

where v = e J/3 — 1, the sum is over the 2'^' spanning subgraphs (V, E') of G, and k(E') is 
the number of connected components (including isolated vertices) of the spanning subgraph 
(V, E'). As a result, the partition function (1.3) is a polynomial in both variables q and v. 
Thus, we can promote (1.3) to the definition of the model, which allows us to consider q as 
an arbitrary complex number. 

We shall here focus on the zero-temperature antiferromagnetic case J = — oo (i.e., v = 
— 1), whose restriction to q G N can be interpreted as a coloring problem. In this case, 
Pg(q) — Zc{q, —1) is known as the chromatic polynomial of G. Its evaluation for a general 
graph is a hard problem. More precisely, Oxley and Welsh have proven that the computation 
of the coefficients of the chromatic polynomial Pg{i) = SL=i a fc<? fc °f a general graph G 
(including bipartite graphs) is #-P hard [4] . 
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In this paper, we will consider only strip graphs of certain regular lattices (i.e., square 
and triangular). For a family of strip graphs {G n } (or, more generally, recursive families 
of graphs [5]), the chromatic polynomial Po n (q) of any member of the family G n can be 
computed from a transfer matrix T and certain boundary condition vectors u and v: 

P Gn = u T -J n -v. (1.4) 

Thus, the computation time grows as a polynomial in n for fixed strip width L. However, 
it still grows exponentially in the strip width (due to the L dependence of dimT). It is 
therefore of interest to devise new and efficient algorithms to be able to handle as large 
widths L as possible. 

Remark. It is important to stress that the Potts spin model has a probabilistic interpreta- 
tion (i.e., has non-negative Boltzmann weights) only when q is a positive integer and v > — 1 
(i.e., positive temperature). The Fortuin-Kasteleyn random-cluster model (1.3), which ex- 
tends the Potts model to non-integer q, has non-negative weights only when q, v > 0. In all 
other cases, the model belongs to the "unphysical" regime (i.e., some weights are negative or 
complex), and some familiar properties of statistical mechanics need not hold. For example, 
even for integer q > 2 and real v in the antiferromagnetic range — 1 < v < 0, where the 
spin representation exists and has non-negative weights, the dominant eigenvalue A* of the 
transfer matrix in the cluster representation need not be simple (i.e., the eigenvector may 
not be unique); and even if simple, it may not play any role in determining Zc n (q, v) because 
the corresponding amplitude may vanish. Both these behaviors are of course impossible for 
any transfer matrix with positive weights, by virtue of the Perron-Frobenius theorem. It 
is important to note that the dominant eigenvalues in the cluster and spin representations 
need not be equal, because the former may have a vanishing amplitude (see below). In any 
case, the Potts model can thus make probabilistic sense in the cluster representation at pa- 
rameter values where it fails to make probabilistic sense in the spin representation, and vice 
versa. It is worth mentioning that for q = 4cos 2 (7r/n) with integer n > 3 and planar graphs 
G, the Potts-model partition function Zc{q,v) admits a third representation, in terms of 
a restricted solid-on-solid (RSOS) model [6, 7] (See also Ref. [8] for an recent study with a 
more comprehensive list of references). 

Even though the zero-temperature limit of the Potts antiferromagnet may not have a 
probabilistic interpretation in the Fortuin-Kasteleyn representation [3] for non-integer q, it 
can nevertheless being studied with the standard tools of Statistical Mechanics and CFT 
with appropriate modifications. In particular, we expect that critical "unphysical regions" 
can be described with non-unitary field theories. (See [9, 10] for a more detailed discussion 
on these points.) 

It is well established that antiferromagnetic models in general, and the chromatic polyno- 
mial in particular, are very sensitive to the choice of boundary conditions. Indeed, different 
choices may lead to quite different thermodynamic limits. In our earlier publications [11-13] 
we studied free and cylindrical boundary conditions: free boundary conditions in the lon- 
gitudinal direction and free or periodic boundary conditions in the transverse direction, 
respectively. More recently [9], we considered cyclic boundary conditions: namely, free 
boundary conditions in the transverse direction, and periodic boundary conditions in the 
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longitudinal direction. In particular, we observed that in the latter case, the phase di- 
agram of the triangular-lattice chromatic polynomial was more involved than for free or 
cylindrical boundary conditions. In this paper we shall deal with fully periodic boundary 
conditions (i.e., periodic in both directions). Periodic boundary conditions in the longitudi- 
nal direction (i.e., toroidal and cyclic) can be easily implemented in the spin representation 
as Pc n = trT™. In the Fortuin-Kasteleyn representation, this implementation is less obvious 
as shown in Ref. [9]. For the triangular-lattice model with periodic boundary conditions on 
the transverse direction we should also use some additional tricks [11, 13, 14]. 

The results existing in the literature about exact chromatic polynomials for strips of 
the square and triangular lattices with toroidal boundary conditions are limited to L < 4. 
The earlier results were obtained by recursive use of the contraction-deletion theorem. These 
works include the pioneering paper by Biggs, Damerell, and Sands [15] (where the solution for 
a square- lattice strip of width L = 2 was presented), and those of Chang and Shrock [16-18]. 
In Ref. [17] the latter authors gave the exact solution for the full Potts-model partition 
function (i.e., with arbitrary v) for square-lattice strips of widths L = 2,3. In Ref. [18] 
the solution for the triangular-lattice chromatic polynomial for width L = 3 was presented, 
while in Ref. [16] they gave the solution for both lattice strips of width L — 4. Finally, in 
2006 Chang and Shrock exhibited a transfer-matrix formalism to compute the full Potts- 
model partition function for regular lattice strips with toroidal boundary conditions [19]. 1 
In particular, they gave the exact partition function for square-lattice strips of widths L < 4, 
and triangular-lattice strips of widths L = 2, 3. 2 In Ref. [19], they also obtained structural 
properties of the partition function and the chromatic polynomial. Additional structural 
properties were discussed by one of us [20] . 

We have several motivations for this work. As explained above, the exact solution of the 
2D q 7^ 2 state Potts model at general v is still unknown, even for the square lattice (and 
a fortiori for more complicated regular lattices), and even in the thermodynamic limit (and 
a fortiori on general finite lattices). Therefore, exact solutions of this model on strips of 
finite fixed width and arbitrary length are valuable in their own right (even in the particular 
case v — — 1). In addition, the chromatic polynomial Pg{q) is an object of great interest to 
mathematicians (see e.g., Ref. [21] for a recent survey). As the computation of the chromatic 
polynomial is #-P hard [4], new and efficient methods for computing this object on large 
graphs are also of great interest to this community. 

The second motivation of this work is to deepen our physical understanding of the critical 
points of the 9- state Potts model in the antiferromagnetic — 1 < v < and unphysical regimes 
v < — 1. Indeed, we can extract this information from our finite- width strips using CFT and 
finite-size scaling (FSS) [22-24]. In particular, we will consider the following way to attain 
the thermodynamic limit. From the chromatic polynomial of a strip of size L x N we 
can compute its chromatic zeros in the complex g-plane. From these zeros it is very hard to 
obtain infinite- volume quantities, as they have strong FSS corrections. Thus, we consider the 
infinite-length limit N — > 00. Then, the chromatic zeros accumulate along certain limiting 

1 Thc transfer-matrix methods of Ref. [19] are essentially the same as ours; but our implementation allows 
us to achieve larger values of the strip width L. 

2 We believe the solution of the chromatic polynomial for triangular-lattice strips of width L = 2 was 
obtained by Chang and Shrock earlier than 2006. In fact, they computed the more involved cases L = 3, 4 
in 2001! 
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curves Bl and around isolated limiting points. This phenomenon is a consequence of the 
Beraha-Kahane- Weiss theorem [25]. Our methods allow us to obtain the exact values of the 
relevant physical quantities in this limit. Their FSS corrections are found to be smaller than 
when both L and N are finite. The extrapolation to the true infinite- volume limit L — > oo, 
using standard FSS techniques, therefore becomes more precise by employing this order of 
limits. 

Remark. Toroidal boundary conditions are rather special compared to the other boundary 
conditions studied in previous papers: namely, free, cylindrical, and cyclic [9,11-13]. Strip 
graphs with the latter boundary conditions are planar] thus, the four-color theorem [26] 
applies. However, strip graphs with the former boundary conditions can be regarded as 
graphs embedded on a torus. In this case, we have only a seven-color theorem. This follows 
from the upper bound obtained by Heawood [27] in 1890 on the chromatic number of a graph 
G embedded on an orientable surface of genus g > 1 



X(G) < H(g) 



I(7 +V /TT48^) 



(1.5) 



combined with the fact that K 7 (with x(^) = 7) can be embedded on a torus, as shown 
in Figure 1. Ringel and Youngs [28] finally proved that the maximum chromatic number 
of any graph G on an orientable surface of genus g > 1 is given by H(g). This graph- 
theoretic difference between graphs embedded on a surface or on a torus may have physical 
implications that are worth studying: in particular, is there any additional structure (e.g., 
phase-transition lines) in the interval q G (4, 7) for the chromatic polynomial of a strip graph 
embedded on a torus? 



The main part of our physical (as opposed to mathematical) understanding of the mech- 
anism for generating chromatic zeros stems from Saleur's analysis of the so-called Berker- 
Kadanoff (BK) phase [29]. This is a massless phase with algebraic decay of correlations. For 
the square lattice, this phase is bounded by the curves 

v±(q) = -2±^~q (1.6) 

and all this region is attracted (at fixed q) by the curve 

v(q) = -y/q, qe [0,4] (1.7) 

It is important to stress four points here: First, the BK phase corresponds to generic values 
of q G (0,4). More precisely, for the Beraha numbers B p 

5 p = 4cos 2 Q, peN, p>2 (1.8) 

the BK phase does not exists due to the vanishing of the amplitude of some of the leading 
eigenvalues and/or cancellations among them [29,30]. Second, the chromatic-polynomial 
subspace v — — 1 intersects the BK phase for q G [0,3]. Thus, because of the attractive 
nature of the BK critical curve (1.7), we can study its conformal properties by considering 
the chromatic-polynomial case in the window q G [0,3]. Third, the upper boundary of the 
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BK phase v+(q) (1.6) can be identified as the antiferromagnetic critical curve [31]. Fourth, 
the exact free energy of the chromatic-polynomial line is not known. However, from our 
previous studies [9,11,12], we can draw a qualitative picture of the chromatic-polynomial 
phase diagram for real q: the system is disordered in the region q G (— oo, 0) U (3, oo), while 
it is critical in the interval q G (0,3). The nature of the system in this critical interval 
depends on the boundary conditions: for free and cylindrical boundary conditions there is 
a single phase [11,12]; but for cyclic boundary conditions, we find two different phases: the 
intervals q G (0, 2), and q G (2, 3) corresponding to two distinct values of a topological order 
parameter (See Section 2.2) [9]. If we define the value q c as the largest value of q such that 
the system is disordered for all q > q c and the system is critical at q = q c [32], then for the 
square lattice we have that q c (sq) = 3 [33]. 

The situation for the triangular lattice is slightly different. This model is solvable on the 
curve [34] 

v 3 + 3v 2 = q (1.9) 

The lower branch of the curve (1.9) can be identified with the lower bound of the BK phase 
V-(q), while the middle branch of (1.9) is expected to play the same role as the curve (1.7) 
[i.e., it is the renormalization-group (RG) attractor governing the BK phase]. Thus, for the 
triangular lattice at least one curve is missing (if we assume a RG flow similar to that of 
the square lattice): the upper bound of the BK phase [i.e., the analogue of the curve v + (q) 
(1.7)], and the antiferromagnetic critical curve. In the limit q — > 0, we have found [10] that 
there is a single phase-transition line separating the (critical) BK phase and the (non-critical) 
high-temperature phase. This line is interpreted as the antiferromagnetic critical curve for 
this model, and its slope is given approximately by 

= -0.1753 ±0.0002 (1.10) 

q=0 

The position of these two curves is an open problem (Some preliminary numerical results 
about the antiferromagnetic critical curve were published in [10, Figure 2]. The full curve 
will be published elsewhere [35]). Theoretical [30] and numerical [10] studies show that the 
BK phase does exist in the triangular-lattice model, and the results suggest that universality 
holds for this phase (i.e., the critical properties of the square- and triangular- lattice BK 
phase coincide). 

Even though some important curves of the phase diagram of the triangular-lattice chro- 
matic polynomial are not known, Baxter [14] found its exact free energy, albeit with peculiar 
boundary conditions that amount to not giving a weight q to clusters with non-trivial homo- 
topy. We comment further on that solution in Sections 5-6 (see also Ref. [13]). The resulting 
phase diagram consists of three different phases in the complex g-plane: the disordered phase 
contained the interval q G (— oo, 0) U (4, oo), and there are two critical phases. One contains 
the interval q G (0, ^(tri)), where ^(tri) is given by 

g 2 (tri) = q B ~ 3.8196717312 (1.11) 

and the third one contains the interval q G ((^(tri), 4). From this phase diagram, we obtain 
that q c (tn) = 4. 

Our previous studies [9,13], gave a similar qualitative picture of that phase diagram. In 
these papers, we have defined qo(L) as the largest value of q at which the limiting curve Bl 
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crosses the real g-axis. We expect that the limit go = hniL^oo qo(L) exists, and it is equal to 
g 2 (tri). 3 In Ref. [13] we reanalyzed Baxter's solution and found that if his three analytical 
expressions for the free energy are correct and exhaustive, the phase transition should not 
occur at ^(tri) = ge, but rather at 



Unfortunately, our numerical data were not conclusive enough to judge which is the cor- 
rect value of the transition point. For free and cylindrical boundary conditions [13], the 
corresponding phase diagram contains only the above three phases. For cyclic boundary 
conditions [9], we found a richer phase structure within the interval q G (0,4). In particular, 
in the interval q G (0,g2(tri)), there should be phase transitions at the Beraha numbers 
q = E>4, B e , . . . < g2(tri) [cf. Eq. (1.8)]. If g2(tri) = ge, then the last such transition occurs 
at q = B u ~ 3.8019377358; but if g2(tri) = g G , then the last transition is at g = B\ 2 = go- 
Finally, the chromatic-polynomial line v — — 1 has a non-void intersection with the BK 
phase for the triangular lattice. In particular, we expect that the intersection should include 
the whole interval q G (0,g2(tri)). Thus, the critical properties extracted from this interval 
are related directly to those of the BK phase. In the interval g G (g 2 (tri), 4), the universality 
class should be different from that of the BK phase. For other boundary conditions, we 
have been unable to extract any meaningful physical information for this interval; so this 
question remains open. In the ferromagnetic regime, it is well-known that results obtained 
from toroidal boundary conditions are closer to the true thermodynamic limit that for the 
other boundary conditions. Although the argument leading to this conclusion may not apply 
for antiferromagnetic or "unphysical" models, it opens an opportunity to try to answer the 
above-posed questions. 

This paper is organized as follows: in Section 2 we show in detail how to obtain the 
chromatic polynomial of a strip with toroidal boundary conditions, and some structural 
properties: dimensionality of the transfer matrix and properties of the amplitudes. Sections 3 
and 4 contain the detailed description (when possible) of the chromatic polynomial for square- 
and triangular-lattice strips, respectively. In these two sections we also describe the limiting 
curve for each regular lattice. The less-mathematically inclined readers can skip most of these 
sections. In Section 5 we discuss the finite-width results and extract the physically interesting 
quantities: phase diagram of the chromatic polynomial for the square and triangular lattices, 
and critical properties of the different phases. Finally, in Section 6, we draw our conclusions. 
In Appendix A we prove two useful combinatorial identities needed in Section 2.2. 

2 Chromatic polynomials with toroidal boundary con- 
ditions 

We want to build the transfer matrix for square- and a triangular-lattice strips of width 
m and arbitrary length n with toroidal boundary conditions. The goal is to write the transfer 
matrix T(m) as a product of two matrices 



g 2 (tri) = g G 



2 + V3 « 3.7320508076 



(1.12) 



T(m) 



H(m) • V(m) 



(2.1) 



3 For the square lattice, on the contrary, we found that go(sq) = <Zc(sq) = 3 [9, 12]. 
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where H (resp. V) represents the horizontal (resp. vertical) bonds of the lattice, and to write 
its chromatic polynomial in the form 

^pxnp(g) = ™ T • H(m) • T(m) n • v id (2.2) 

for some left and right vectors u and v id . The subindex P denotes the periodic boundary 
condition in the corresponding direction. 

2.1 The transfer- matrix formalism for toroidal boundary condi- 
tions 

The main problem in writing a transfer matrix in the Fortuin-Kasteleyn representation 
is the non-local factor q k ( E ') in (1.3) [11]. For free boundary conditions in the longitudinal 
direction (i.e., free and cylindrical), this can be efficiently handled by keeping track of the 
connectivity state v-p of the top layer. This connectivity is related to a certain partition V 
of the single-layer vertex set {1, 2, ... , m} [11, 13]. 

For periodic boundary conditions in the longitudinal direction (i.e., cyclic and toroidal) we 
must keep track of the connectivity state of both the top and bottom rows: each connectivity 
state v-p is now given by a certain partition V of the vertex set of the top row {1,2,..., m) 
and the vertex set of the bottom row {1', 2', . . . , m'}. In Figure 2 we show an example of a 
square-lattice strip of size 4x2. As for cyclic boundary conditions [9], initially the top and 
bottom rows are identical. Then, we enlarge the strip by adding new layers via the action of 
the transfer-matrix operator T(m) exclusively on the top-row vertices (while the bottom-row 
sites remain unchanged). At the end, when we have built a strip of n + 1 rows, we identify 
the top and bottom rows. 

Let us now explain in more detail how to obtain the chromatic polynomial for a square- 
lattice strip of width m, length n, and toroidal boundary conditions. 

We characterize the state of the top and bottom rows by a connectivity state v-p, which 



is associated to a partition V = {V±, V 2 , • • • , Vk] of the set {1,2,..., m, 1', 2', . . . , m'}. We 
then introduce the operators 

P x = vI + D x (2.3a) 

Q X ,y = I + VJ x ,y (2.3b) 



acting on the connectivity space {vp}. Here I is the identity, is the detach operator 
that detaches site x from the block it currently belongs to, and S x>y is the join operator that 
amalgamates the blocks containing sites x and y, if they were not already in the same block. 
(Further details can be found in Ref. [11].) 

The matrices V and H are defined in terms of the above operators as follows: 

(m— 1 \ 
J] Qm+iJ • Qm,i (2.4a) 

m 

V sq (m) = ( 2 - 4b ) 
i=i 

Notice that periodic boundary conditions have been implemented by adding the operator 
Q mi i in (2.4a). We can similarly define a matrix \-\' sq acting only on the bottom row. 
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Then, the chromatic polynomial for a square-lattice strip of size m x n with toroidal 
boundary conditions can be written as in (2.1)/(2.2), where the two vectors are those of 
cyclic boundary conditions [9]. The vector v id denotes the partition 

v id = {{l,l'},{2,2'},...,{m,m'}} (2.5) 

(i.e., we start with the top and bottom rows identified). The left vector u T acts on a 
connectivity state v v as 

u T -v v = q\ v '\ (2.6) 
where \V'\ denotes the number of blocks in the connectivity state v' v obtained from v v by 

v' v = (ilvj - vv (2-7) 

Thus, u T acts on v v by identifying the top and bottom rows, and then by assigning a factor 
of q to each block in the resulting partition. 

The definition of the transfer matrix for a triangular-lattice strip of width m, length 
n, and toroidal boundary conditions is more involved because of the periodic boundary 
conditions in the transverse direction. As explained in Refs. [11,13], in this case one has to 
consider a triangular lattice of width m + 1 with free boundary conditions in the transverse 
direction, and then identify the columns 1 and m + 1. (See Figure 3 for an example of size 
4x2.) Then, the matrices H and V take the form 

Htri(m) = Q M+1 j • J m+1 ,i (2.8a) 

Vtri(m) = (f]P,-Q M+1 ] -P m+1 (2.8b) 



,?=i 



where the operators Q^+i in (2.8b) represent the oblique bonds e = + 1), (i + 1, j)). 
With this definitions, the chromatic polynomial can be written as in (2.1)/(2.2) with the 
same right and left vectors as for the square-lattice case (2.5)/(2.6). 

A further simplification can be obtained if we take into account the fact that v = — 1. 
Then, H is a projector (i.e., H 2 = H), and we can use instead of the transfer matrix T and 
the basis vectors v v , the modified transfer matrix T' 

l\m) = H(m) • V(m) • H(m) (2.9) 

and the modified basis vectors 

w v = H(m) • v v (2.10) 

Note that w v — if V has any pair of nearest-neighbor sites of the top row in the same 
block. Thus, the space {w-p} has lower dimension than the original connectivity space {v-p}. 

Notice that we start with the state = H • v^, and that at the end we identify the top 
and bottom rows. Then, it is very useful to work directly with the modified connectivity 
basis 

w v = H'(m) • H(m) ■ v v (2.11) 
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This choice implies that w v = also if V contains any pair of nearest-neighbor sites of 
the bottom row in the same block. For simplicity, we will drop hereafter the prime in the 
modified transfer matrix (2.9) and the hat in the modified basis vectors (2.11). Then, the 
chromatic polynomial can be written as 

Pra^nM) = U T -J(m) n -W id (2.12) 

Finally, to lighten the notation, we will write the connectivity states v v using Kronecker 
delta-functions. For instance, V{{i,2},{i',3'},{3},{2'}} will be written as Si^Sy^i. The identity 
connectivity state (2.5) will be denoted simply as 1. 

2.2 Structural properties of the chromatic polynomial 

Our first goal is to compute the number C m ^ m of non-crossing connectivities in the interior 
of an annulus with m points on the inner rim and m points on the outer rim (See Figure 4). 
This problem can be reduced to a simpler one by introducing the quantities d(m, £). Given m 
points on a circle and a non-crossing connectivity V exterior to that circle, we can compute 
the number of ways d(m, £, V) we can connect that particular connectivity V to infinity by 
means of £ mutually non-connected paths (bridges). Then, the d(m,£) are the sum over all 
possible connectivities V 

d(m,e) = ^2d(m,e,V) (2.13) 
v 

The values of d(m,£) are given by the following formula [19] 4 

f C m ii£ = 

if 1= 1 

(2.14) 

if 2 < £ < m 
if £ > m 

where C m is the Catalan number 

C„ = ^- M (2.15) 
m+l \m j 

which gives the number of non-crossing connectivities with m points on a circle (or on a 
line) . 

The number of connectivities C m ^ m can be obtained easily from the d(m,£) (2.14) by 
using the £ paths as bridges and by matching two such bridge connectivities with the same 
value of < £ < m 

m 

C m , m = C 2 m + Y,td(m,£) 2 (2.16) 
l=\ 

4 The numbers d(m,l) (2.14) are denoted as n z ,tor(m, i) by Chang and Shrock [19, Eqs. (2.19)-(2.22)], 
and as n tor {m,l) in Ref. [20, Eq. (2.1)]. 



d(m, 



2m- 1 
m — 1 

2m 
m — £ 
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If we insert (2.14) for the d(m,£), then we obtain a closed expression for C m) m- 

C m „= ro '~^ t 1 +B r aro ~ 1 V P.") 
m + 1 \ m J 

where we have used the identity (A. la) proven in Appendix A. The asymptotic growth of 
C m ,m as m — > oo is exponentially fast with a large constant: 



c 16 



/ m,m ~ 



An 



I - ( - 

m 



(2.18) 



In practice, we build the C m ^ m non-crossing connectivities in a recursive way. We start 
with all the C m _i ;m _i non-crossing connectivities on an annulus with m — 1 points on its inner 
and outer rims, and then we add a new site on the inner rim. The number of non-crossing 
connectivities for this annulus is obtained again by matching the ^-bridge connectivities 



m—1 



C m , m _! = C m C m - 1 + Y t ^d(m,£)d(m-l,e) (2.19) 

t=i 

Then, we add a new site on the outer rim, and we obtain the C m ^ m connectivities (2.17). A 
closed formula for C m>m _i (2.19) can be obtained if we use the combinatorial identity (A. lb) 
proven in Appendix A. After some algebra, we obtain 

2m 4 — 6m 3 + 11m 2 + 3m — 4 / 2m — 2 X 2 



C m,m-1 = „ ' ' 1 (2.20) 

2m 2 (m+l) \m — l J 

valid for all m > 2. This formula also grows exponentially fast in m in the large-m limit: 



c 16 



16tt 



1 + 



as m — > oo (2-21) 



In Table 1 we show the values of C m ^ m and C m) m-\ for 1 < m < 7. 

The condition v = — 1 implies that connectivities V containing nearest-neighbor sites in 
any sub-block do not contribute. The number of non-crossing non-nearest-neighbor connec- 
tivities for an annulus with m points on its inner and outer rims will be denoted by D m m . 
This number gives therefore the dimension of the modified transfer matrix (2.9) for a square- 
or triangular-lattice strip of width m and toroidal boundary conditions. These numbers are 
displayed for 1 < m < 7 in Table 1. 

As in previous works [9, 11-13], we can reduce this dimension by considering the equiv- 
alence classes modulo symmetries of non-crossing non-nearest-neighbor connectivity states. 
For the triangular (resp. square) lattice, we should consider equivalence classes modulo rota- 
tions (resp. rotations and reflections with respect to any diameter of the strip). The number 
of generated equivalence classes modulo rotations (resp. reflections and rotations) of non- 
crossing non-nearest-neighbor connectivity states for a strip of width m will be denoted by 
TriTorus'(m) (resp. SqTorus'(m)). The adjective "generated" means that we only take into 
account those classes of connectivity states that can be produced by applying the various 
operators P, and Qjj on the right vector io id . This distinction is relevant for the square 
lattice, as there are perfectly legal non-crossing non-nearest-neighbor partitions that cannot 
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be generated when we start from the right vector io id . (See e.g., the discussion on this is- 
sue in Ref. [20, end of Section 2.1].) The simplest example corresponds to to = 2 and the 
connectivity state 52,i"5i,2'- In general, out of the m possible connectivity states with £ = m 
bridges, only the state v id (2.5) can be realized on a square-lattice strip. The other to — 1 
non-crossing connectivity states cannot be realized because of the lattice structure. However, 
on the triangular lattice, which has a larger set of "vertical" edges, one can produce those 
m — 1 connectivity states without violating the non-crossing condition. 

Finally, as for cyclic boundary conditions, for a given width to we find many repeated 
eigenvalues. Thus, we define TriTorus(m) (resp. SqTorus(m)) as the number of distinct 
eigenvalues (with non-zero amplitude) for a triangular-lattice (resp. square-lattice) strip of 
width m with fully periodic boundary conditions. In Table 1 we have displayed the quantities 
TriTorus', SqTorus', TriTorus and SqTorus for to < 7 (when available). 

The symbolic computation of the transfer matrix and the ia t and vector has been 
carried out using a mathematica script for the smallest widths 2 < to < 4, and a C program 
for the larger widths to < 7. Indeed, for to < 4, both methods perfectly agree. We refer to 
Refs. [9, 12, 13] for details about the symbolic computation. For larger widths to < 15, we 
have computed numerically the transfer matrix using another C program. 

Remarks. 1) With our transfer- matrix approach we have checked all previous results by 
Biggs, Damerell and Sands [15], and Chang and Shrock [16-19]. 

2) Even though the final dimensions of the transfer matrices are TriTorus' (to) or SqTorus' (m), 
the actual computations are carried out in the larger space of non-crossing connectivities (of 
dimension C m>m ). This fact limits the maximum width we are able to handle. 

3) We have cross-checked our results using the trivial identity 

(?) (?) (2-22) 

for all 2 < to, n < 7. 



2.3 Amplitudes for toroidal boundary conditions 



The transfer matrix for a strip with toroidal boundary conditions has a block structure. 
This property also holds for cyclic boundary conditions [9], and in both cases the origin is 
the same. Given a certain connectivity V of the top and bottom rows (see e.g., Figure 4), 
we can regard it as a bottom- row connectivity Pbot, a top- row connectivity P t op joined by I 
bridges. Notice that each top-row block can be connected to a bottom-row block by at most 
one bridge, and vice versa. Thus, the number of bridges is an integer < £ < to. 

As explained in Ref. [9] , the full transfer matrix has a lower-triangular block form (if we 
order the connectivity states in decreasing order of bridges). 



T(to) 



/ T 

T m —l m 







T 

J- i 



\ 



0,m 



m—l,m—l 



-0,m-l 



\ 




Tr 



(2.23) 



o,o / 



The reason is the application of the H and V operators cannot increase the number of 
bridges of a given connectivity state. Furthermore, as the transfer matrix does not act on 
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the bottom-row connectivity and they act on the space of non-crossing connectivities, each 
diagonal block Tgj has also a diagonal-block form 



/ Tt 



Ti,t,i 



\ 







(2.24) 



V 




where each sub-block Tgjj is characterized by a certain bottom-row connectivity Vbot and 
a position of the £ bridges. Its dimension is given by the number of top-row connectivities 
that are compatible with Vbot, and the number £ and relative positions of the bridges. The 
number of blocks Ng (0 < £ < m) for square- and triangular-lattice strips of width m are 
displayed in Tables 2 and 3, respectively. 

In particular, as for cyclic boundary conditions [9], this structure means that the char- 
acteristic polynomial of the full transfer matrix T can be factorized as follows 



In practice, for each strip width 2 < m < 7, we have performed the following procedure: 

1. Using a C program, we compute the full transfer matrix T(m) and the left u and 
right vectors W\&. These objects are obtained in the basis of non-crossing non-nearest- 
neighbor classes of connectivities that are invariant under the right symmetries (de- 
pending on the lattice structure). The dimension of this transfer matrix is SqTorus'(m) 
or TriTorus'(m). With these objects, we can compute the chromatic polynomial of any 
strip of finite length n using Eq. (2.12). 

2. Given the full transfer matrix T(m), we compute the the SqTorus(m) or TriTorus(m) 
distinct eigenvalues Xi(q). To do this, we split the transfer matrix into blocks Tg^j, 
each block characterized by a bottom-row connectivity and the number £ and position 
of the bridges (modulo the corresponding symmetries). By diagonalizing these blocks 
we obtain the whole set of distinct eigenvalues \i(q) and their multiplicities ki. 

3. The final goal is to express the chromatic polynomial in the form 



where the {Aj} are the distinct eigenvalues of the transfer matrix T(m) and the {a{\ 
are some amplitudes we have to determine. The calculation of the amplitudes can 
be achieved by solving SqTorus(m) or TriTorus(m) linear equations of the type (2.26), 
where the l.h.s. (i.e., the true chromatic polynomials) have been obtained via Eq. (2.12) 
(i.e., by using the full transfer matrix and vectors). 5 

5 In practice, to obtain the amplitudes at(q) we assume that they are polynomials in q. The order of the 
polynomial cti(q) corresponds to the largest £ sector the eigenvalue Aj belongs to. Unlike for cyclic boundary 



rn 



x(v = n ii^tw 



(2.25) 




(2.26) 
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We have used this procedure to compute the transfer matrix up to width m = 7 for both 
the square and the triangular lattices. The next cases were unmanageable with our current 
computer resources. 6 The number of distinct eigenvalues found in each block can be 
found in Tables 4 and 5 for the square and triangular lattices, respectively. 

In Ref. [9], we found several properties of the eigenvalues Aj and their corresponding 
amplitudes ai for cyclic boundary conditions. In particular, that all eigenvalues within a 
block T^i have the same amplitude (see below), and that any two distinct blocks T^g 
and T(t£i (with £ ^ £') have no common eigenvalues. These properties were essentially due 
to the presence of a quantum-group symmetry in the Potts model with cyclic boundary con- 
ditions. These properties were also used to simplify the algorithm to compute symbolically 
the transfer-matrix eigenvalues Aj. Unfortunately, this quantum-group symmetry is broken 
when toroidal boundary conditions are considered, and none of these properties hold in our 
case. 

Chang and Shrock [19] have found that the chromatic polynomial for toroidal boundary 
conditions (2.26) has a structure similar (but not identical) to that for cyclic boundary 
conditions [9,29,36,37]. For a square- or triangular-lattice strip of width m, length n and 
toroidal boundary conditions, the chromatic polynomial can be written as 

P mpX n P (g) = EX>f A ^( m ) n ( 2 - 27 ) 

e=o j=i 

where the amplitudes b*p are polynomials in q of order £, and in contrast with cyclic boundary 
conditions, are not the same within a given £ sector. The eigenvalues A^j(m) are algebraic 
functions of q and, in contrast with cyclic boundary conditions, a given eigenvalue can appear 
in more than one £ sector. 

Chang and Shrock [19] found a family of basic amplitudes given by 

6 (o) = a (o) = x 

= = q -l 

and for £ > 2, by 

&W = a W_ a H) + (-i/ Q W (2.29a) 
= <f + B-l) 5 ^~ ( 2i ~ j ) ^ + (9 " !) ( 2 - 29b ) 



conditions, a given eigenvalue for toroidal boundary conditions can appear in several sectors (each of them 
characterized by a different value of I. See text.) Once the Ansatz for the amplitudes is fixed, we solve the 
linear equations for the numerical coefficients in a^q) using as many chromatic polynomials P mp xn P (q) as 
needed. After finding the solution, we check that it can reproduce chromatic polynomials not used in the 
determination of the amplitudes. 

6 Thc largest transfer matrix we have computed corresponds to the m = 7 triangular lattice of dimension 
40996. We needed ~ 4.3Gb of RAM memory and ~ 140 hours of CPU to perform this computation. In 
addition, the raw mathematica output file was very large ~ 355Mb. 
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(2.28a) 
(2.28b) 



where the coefficients are the amplitudes obtained for cyclic boundary conditions [29] 



a 



2 



3=0 



(2.30a) 
(2.30b) 



where C/ n (x) is the Chebyshev polynomial of second kind. The first non-trivial are [19, 
Eqs. (2.1)-(2.4)]: 





= g 2 


-3g + l 




(2.31a) 


b m 


= g 3 


- 6q 2 + 8q - 1 




(2.31b) 


b w 


= g 4 


- 8g 3 + 20g 2 - 15g + 1 




(2.31c) 


6 (5) 


= g 5 


- 10g 4 + 35g 3 - 50g 2 + 24g - 1 




(2.31d) 


6 (6) 


= g 6 


- 12g 5 + 54g 4 - 112g 3 + 105g 2 - 


35g + 1 


(2.31e) 


6 (7) 


= g 7 


- 14g 6 + 77q 5 - 210g 4 + 294g 3 - 


196g 2 + 48g - 1 


(2.31f) 



Finally, Chang and Shrock argue that the coefficients b^p cannot be equal to the above b^. 
For m = 2, 3 they find that the coefficients b^ and b^ split into two coefficients (up to a 
positive integer factor): 7 





b? = \q(q-3) 


(2.32a) 




= i(g-l)(g-2) 


(2.32b) 




b? = i(g-i)(g 2 -5g + 3) 


(2.32c) 




bf = \q{q-2){q-A) 


(2.32d) 


In our numerical study, 
amplitudes are given by 


we have found that this also true for m = 4,5. 


The basic 


&! 4) 


= ^g(g-2)(g-3) 2 


(2.33a) 




= ^g(g-l)(g-3)(g-4) 


(2.33b) 


6 (4) 


= I( g -i)( g 3_ 7g 2 + 14g _ 4) 


(2.33c) 


b? 


= i(g 5 -10g 4 + 35g 3 -50g 2 + 24g-5) 
5 


(2.33d) 




= ig(g-l)(g-2)(g-3)(g-4) 




(2.33e) 



Note that the expression for 4 3) (2.32d) is twice the expression given by Chang and Shrock [19, 
Eq. (2.36)]. 
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In general, an eigenvalue Aj which only appear in a single block has an amplitude given 
by a positive integer multiple of either (2.28)/(2.31) or (2.32)/(2.33). All the amplitudes 
and b^ (2.31) are simply polynomials in q of order I. When an eigenvalue \ appears in 
more than one block T^, then its amplitude is just a certain linear combination (with integer 
coefficients) of the corresponding amplitudes given above. Thus, the order of this polynomial 
amplitude is just the largest value of £ involved. 

Remarks. 1. After we empirically obtained the amplitudes (2.33), one of us obtained the 
general formula for the different amplitudes b^p [20]. Indeed, our results fully agree with the 

general formula. The number of distinct amplitudes b^ for a given £ was shown to be the 
number of integer divisors of £ [20]. 

2. As we shall discuss in the next Section, the sub-block T LjL of the square-lattice strip 
of width L contains the eigenvalue A = (— 1) L with amplitude b^ (2.31). Indeed, each b^ 
can be written as a linear combination of the amplitudes b^ (2.32)/(2.33): 

= & j 2 )+4 2 ) (2.34a) 

6(3) = bf ] +2bf (2.34b) 

6^ = b { ? +2b { f> + bf (2.34c) 

6 (5) = b f +Ah f ( 2 .34d) 

This fact was first noted by Chang and Shrock [19]; and has been proven in Ref. [20, c.f., 
Eq. (4.7)] by showing how many times each b, appears in a given 

3. In our previous study of the chromatic polynomial for cyclic boundary conditions [9], 
we could only obtain the exact chromatic polynomial up to widths m — 6. The reason was 
that the computation of the characteristic polynomial of a symbolic matrix of dimension 
> 100 was unfeasible with mathematica's standard built-in functions (based essentially on 
Gaussian elimination). In this paper, we have been able to overcome this problem by using a 
different division-free algorithm, which computes the characteristic polynomial of a matrix 
with elements belonging to a ring (not to a field). We have used a mathematica implemen- 
tation [38] of the Samuelson-Berkowitz-Abdeljaoued algorithm [39]. 8 With this algorithm 
we have been able to deal with symbolic matrices of dimensions up to 308. Unfortunately, 
for L = 7 we need to compute the symbolic characteristic polynomial of much larger ma- 
trices (of dimensions 483 and 532). These computations are beyond our current computer 
facilities, even with the help of the above-described improved algorithm. 9 Further details 
will be published elsewhere. 



3 Square- lattice numerical results 

In this section we will analyze the results for the square- lattice strips of widths 2 < L < 7. 
We have checked our results using (in addition to the already known cases L < 4 [16,17,19]) 
the identity (2.22). 

8 We thank Alan Sokal for bringing to our attention several important references on this topic. 

9 Thc computation of the characteristic polynomial of the sub-block of dimension 308 took approximately 
two months of CPU time. 
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For each L, we give a detailed (when possible) description of the eigenvalues A; and 
amplitudes ctj. We have included the already known cases L < 4 for completeness and, 
more importantly, to show in detail how our method works. In particular, for L = 2, 3 we 
give the explicit expression of the transfer matrix and the left and right vectors. The exact 
expressions for the new cases L = 5, 6 are very lengthy, so we refrain from writing down all 
the needed formulae. Instead, we have included a mathematica file named transf er6_sqT.m 
available as part of the on-line version of this paper in the cond-mat archive at arXiv . org. 
For L = 7, even though we have obtained the symbolic form of the transfer matrix and the 
left and right vectors, we have been unable to compute some of the eigenvalues (because of 
the large dimensionality of some sub-blocks), and thus, the expression of the amplitudes. 

Finally, for each strip, we compute the corresponding limiting curve and isolated limiting 
points (where the chromatic zeros accumulate in the infinite-length limit), and analyze their 
main properties. 



3.1 L = 2 



The chromatic polynomial for this strip is exactly the same as the chromatic polynomial 
for a square- lattice strip of with L = 2 with cyclic boundary conditions [9]: in the basis 
{^1,1^2,2', <5i,i' + $2,2', 8i,2> + $2,i', 1}, h takes the form 



T(2 T ) = 



/ 


i 





\ 


-1 




2-q 1 
1 2-q 






V 


i 


2(q-2) 2(q-2) 


q 2 


-3q + 3 J 



(3.1) 



where we show by vertical and horizontal lines the lower-triangular block structure of this 
matrix. The right w id and left u vectors are given by 

u T = (A,2A,0,A), 



wf d = (1,0,0,0) 



As in [9] , we find that the three diagonal blocks for I 
eigenvalues and amplitudes: 

• £ = 0: The eigenvalue is A = q 2 — 3q + 3 with amplitude 



with A — q(q— 1) (3.2a) 

(3.2b) 

- 0, 1, 2 bridges give the following 



1. 



£ = 1: There are two eigenvalues Ai ; i 
amplitude = q — 1. 



1 — q and Ai j2 = 3 — q, with the common 
£ = 2: The eigenvalue is A 2 = 1 with amplitude = q 2 — 3q + 1. 



Each diagonal block is characterized by the bottom-row connectivity Vbot = 1- 

In Figure 5(a) we have shown the chromatic zeros for the strips of lengths n — 10 and 
n = 20. We also show the limiting curve £>2 where the chromatic zeros accumulate in the 
infinite-length limit. 

The limiting curve £> 2 splits the complex g-plane into four regions. The regions inside the 
outer parts of £> 2 are characterized by a dominant eigenvalue belonging to the £ = 1 sector, 
while in the outer region the eigenvalue Ao dominates. 

The curve £> 2 crosses the real q axis at two points: q = and q — 2. This latter point is 
in fact a multiple point. Finally, there is a single isolated limiting point at q — 1. 
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3.2 L = 3 



The chromatic polynomial for this strip was computed by Chang and Shrock [17] (in fact, 
they computed the full Potts-model partition function [19]). 

The full transfer matrix T(3t) has dimension eight. In the basis {^i,i'52,2'^3,3', ^i,2'^2,i' + 
• • • , ^i,2'^3,3' + • • • , 5i 5 i/52,2' + • • • , 5i,2'5i,2' + . . . , #i ; 2' + • • • , ^i,i' + • • • , 1} (where the dots ". . ." 
mean all possible states equivalent under reflections and/or rotations), it takes the form 



T(3 T ) 



/ 


-1 




















\ 







S 3 








-2 



















S3 


-1 


-1 













i 





-2 


S 3 



















-1 


-1 





S3 
















-s 4 


-s 4 


1 


7-2q 


S 3 - S 


S3 







-1 





-2s 4 


-2s 3 


2 


2s 3 


-s 





V 


1 


3s 3 


6s 3 


3s 3 


6s 3 


6S 


3S 


Ao / 



(3.3) 



where we by vertical and horizontal lines show the lower-triangular block structure of this 
matrix, and 



Sfc 

S 
A 



q — k 
2 5q + 7 



= q 



= q 3 - 6q 2 + Uq - 13 



The vectors u>id and u take the form 

■ T -= (A,0,0,3A,0,0,3A,A), 



u 



with A — q(q— l)(q - 2) 



w 



T 



(3.4a) 
(3.4b) 
(3.4c) 

(3.5a) 
(3.5b) 



M(1 - (1,0,0,0,0,0,0,0) 

The eigenvalues can be obtained from the diagonal blocks Tg 7 g. Each block Tgj is char- 
acterized by the trivial bottom- row connectivity state Vhot = 1- We find the following 
eigenvalue structure for T(3t) (3.3): 

• £ = 0: There is a single eigenvalue Aq (3.4c) with amplitude = 1. 



£ = 1: We find two eigenvalues: Ai 5 i = — (q 2 — 7q + 13) with amplitude 
and Ai >2 = -(q - 2) 2 with amplitude 2b {1) = 2{q - 1). 



1 



£ = 2: The four eigenvalues are A2,i — q — 1, \2,2 — Q~ 2, A2,3 
The amplitudes are given respectively by b^\ 2b^ 



q — 4, and A 2 ,4 = g — 5. 



2b 2 *>, and bf ] [c.f., (2.32)]. 



• £ = 3: The eigenvalue is A 3 = —1 with amplitude b^ (2.31b). 

In Figure 5(b) we have shown the chromatic zeros for the strips of lengths n — 15 and 
n = 30, as well as the limiting curve £> 3 where the chromatic zeros accumulate in the infinite- 
length limit. This curve splits the complex g-plane into three regions. The region containing 
the real segment (0, 2) is dominated by an eigenvalue from the £ = 1 sector. The region 
containing the real segment (2, 3) is dominated by the £ = 2 sector. The rest of the complex 
g-plane is dominated by the £ = eigenvalue Ao (3.4c). 

The curve £> 3 crosses the real q axis at three points: q = 0, q = 2, and q = 3. We also 
find a pair of complex conjugate T points at q ~ 2.4692972375 ± 1.4013927796 i. Finally, 
there is a single isolated limiting point at q — 1. 
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3.3 L = 4 



The chromatic polynomial for this strip was computed by Chang and Shrock [16]. The 
full transfer matrix T(4t) has dimension 68. In this case, we restrict ourselves to report the 
eigenvalue structure: 

• t = 0: This block has dimension five with two sub-blocks of dimensions three and 
two. There are three distinct eigenvalues: Ao,i = (q — l){q — 3) and the solutions of a 
second-order equation [16, Eq. (2.12)]. All of them have the amplitude = 1. 

• t = 1: This block has dimension 22, and it contains three sub-blocks of dimensions 
six and eight (two of them). We find eight different eigenvalues. Two of them are 
the solutions of the second-order equation [16, Eq. (2.13)] with amplitudes 2b^ = 
2(q — 1); The other six eigenvalues come from two distinct third-order equations [16, 
Eqs. (2.18)/(2.19)], with amplitudes =q-l. 

• £ = 2: This block has dimension 33 containing three sub-blocks of dimensions 8, 12, 
and 13. Among these 33 eigenvalues, we find 15 distinct eigenvalues. Three of them are 
simple: A 2 ,i = q 2 — 5g + 5 with amplitude A 2j2 = q 2 — 5q + 7 with amplitude bf \ and 
the special one A( 2i 3) = 2— q which also appears in the £ = 3 sector (see below). The next 
six eigenvalues are solutions of three second-order equations [16, Eqs. (2.14)-(2.16)], 
with amplitudes bf \ 2bf\ and bf \ respectively [c.f., (2.32)]. The last six eigenvalues 
come from two third-order equations [16, Eqs. (2.20)/(2.21)], with amplitudes bf^ and 
2bf\ respectively. 

• I = 2,3: The eigenvalue A( 2;3 ) = 2 — q appears in these two blocks with an amplitude 
(q - 2)(4g 2 - 13q - 3)/6 = bf + 2bf . 

• t = 3: This block has dimension seven: in addition to A( 2j 3), we find six other eigen- 
values: A3 ; i = 5 — q with amplitude b\ [16, Eq. (2.8)]; A3 i2 = 4 — q with amplitude 
2bf [16, Eq. (2.7)]; A 3i3 = 3 - q with amplitude 2bf ] [16, Eq. (2.6)]; A 3i4 = 1 - q 
with amplitude bf 1 ^ [16, Eq. (2.4)]; and \3,5,± = 3 — q ± with common amplitude 
2bf [16, Eq. (2.17)]. 

• I = 4: This one-dimensional block gives A4 = 1 with amplitude b^ (2.31c). 

Thus, for this strip we find 33 different eigenvalues. All of them but A( 2;3 ) = 2 — q, belong 
to a single £ sector. 

In Figure 5(c) we have shown the chromatic zeros for the strips of lengths n = 20 and 
n = 40, as well as the limiting curve B4. This curve splits the complex g-plane into five 
regions. The region containing the real segment (0, 2) is dominated by an eigenvalue from 
the i = 1 sector. The region containing the real segment (2, <?o(4)) with qo(4) m 2.7827657401 
is dominated by the £ = 2 sector. The rest of the complex g-plane (i.e., the other three 
regions) is dominated by the £ = sector. Notice that the special eigenvalue A( 2)3 ) = q — 2 
is not dominant in any open set of the complex g-plane. 

The curve B4 crosses the real q axis at three points: q = 0, q = 2, and q = (?o(4) ~ 
2.7827657401. We also find three pairs of complex conjugate T points at q ~ 2.4184666383 ± 
1.7281014476 i, q w 2.6228997762 ± 1.5279208067 i, q m 2.7784542034 ± 0.9708262424 i. 
Finally, there is a single isolated limiting point at q — 1. 
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3.4 L = 5 

The full transfer matrix T(5t) has dimension 347. The relevant eigenvalues can be 
extracted from the diagonal blocks T«(5t). We find the following eigenvalue structure: 

• £ = 0: This block has dimension six. In the basis {Sr t y + . . . , 5i',3'5i,3 + . . . , 5i' ,3^1,4 + 

. . " mean all equivalent states under 

\ 



. . . , ^i',3'i52,4 + . . • , 1, 5i,3 + . . .} (where the dots 
rotations and reflections), it takes the form: 



2o,o(5i 
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5T 2 





(3.6) 



where s k and 5 are defined in (3.4a)/(3.4b), and 



T, = 
T 2 = 

n = 



g 5 - 10g 4 + 45g 3 - 115g 2 + 169g - 116 



Q 



9g 3 + 34q 2 - Q3q + 47 



<f - 9q 2 + 29q + 32 



(3.7a) 
(3.7b) 
(3.7c) 



This matrix is block diagonal: the first four- dimensional block corresponds to the 
bottom-row connectivity Pbot = ^i',3', while the second two-dimensional block corre- 
sponds to the trivial bottom- row connectivity V^ot — 1- 

We find four distinct eigenvalues coming from two second-order equations: 

= ^ 2 - £(2q 3 - 13q 2 + 28q - 19) + g 6 - 13g 5 + 69g 4 - 191g 3 

+292g 2 - 236g + 79 (3.8a) 

= £ 2 - £(g 5 - 10g 4 + 46g 3 - 124g 2 + 198g - 148) + q 8 - I9q 7 + 159g 6 

-767g 5 + 2339g 4 - 4627g 3 + 5800g 2 - 4212g + 1362 (3.8b) 

Their amplitudes are 26 (0) = 2 and 6 (0) = 1, respectively. 

• £ = 1: This block has dimension 67 and it contains four sub-blocks of dimension 14 
(three of them) and 25. We find 15 different eigenvalues: An^) = — {q 2 — 4g + 3) also 
appears in the i = 2 block (see below), and the rest come from polynomial equations 
of order ten (with amplitude 2b^ = 2(q — 1)) and four (with amplitude = q — 1). 

• I = 1,2: The eigenvalue A(i ;2 ) = — (q 2 — 4g + 3), appears in these two blocks with 
amplitude (q + l)(q - 2)/2 = + bf\ 

• t = 2: This block has dimension 190 and it contains five sub-blocks of dimensions 
31, 33 (three of them), and 60. We find 36 distinct eigenvalues, including A(i i2 ). Two 
eigenvalues come from the second-order equation: 



i 2 + i{2q 2 - \\q + 16) + g 4 - llg 3 + 44g 2 - 77g + 51 







(3.9) 
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with amplitude b 2 (2.32b). The other eigenvalues are the solutions of polynomial 
equations of order four, five, and twelve (two of them). Their amplitudes are are 
&S 2) , 2b ( ?\ and 2b ( 2 2 \ respectively [c.f., (2.32)]. 

• £ = 3: This block has dimension 72, and it contains three sub-blocks of dimension 
24. The 24 distinct eigenvalues come from polynomial equations of order three (two of 
them), six, and twelve. The corresponding amplitudes are 2b 2 °', b?', 2bf, and 2b y 2 >, 
respectively [c.f., (2.32)]. 

• t = 4: This block has dimension 11. Three eigenvalues are rather simple: A^i = q — 5 
with amplitude b^f 1 ; A^ = q — 3 with amplitude 2b^ ; and A^ = q — 1 with amplitude 
b^ [c.f., (2.33)].. The next four are solutions of two second-order equations 

= £ 2 _£(2g-5) + g 2 -5g + 5 (3.10a) 
= e-a^-7) + q 2 -7q + ll (3.10b) 

with amplitudes 2b^ and 2b^ , respectively. The last eigenvalues come from a fourth- 
order equation with amplitudes 2b^ ■ 

• £ = 5: This one-dimensional block gives A 5 = — 1 with amplitude b^ (2.31d). 

For this strip we find 90 different eigenvalues. All of them but one belong to a single £ 
sector. The only exception is the eigenvalue A(i ;2 ) = — (g 2 — 4g + 3) which appears in the 
sectors £ = 1,2. 

In Figure 5(d) we have shown the chromatic zeros for the strips of lengths n = 25 and 
n = 50, and the limiting curve 3%. This curve splits the complex g-plane into three regions. 
The region containing the real segment (0, 2) is dominated by an eigenvalue from the £ = 1 
sector. The region containing the real segment (2, 3) is dominated by the £ = 2 sector. 
The rest of the complex g-plane is dominated by the £ = sector. Notice that the special 
eigenvalue An ; 2) is not dominant in any open set of the complex g-plane. 

The curve £>5 crosses the real q axis at three points: q = 0, q = 2, and q = 3. We 
also find two pairs of complex conjugate T points at q ~ 2.2834116867 ± 2.0404211832 i, 
q Ki 2.4349888915 ± 1.9935503628 i. We also find a pair of complex-conjugate endpoints at 
q ss 2.5034648020 ± 2.0851731763 i. Finally, there is a single isolated limiting point at q = 1. 

3.5 L = 6 

The full transfer matrix T(6t) has dimension 2789. The relevant eigenvalues can be 
extracted from the diagonal blocks T«(6t)- We find the following eigenvalue structure: 

• £ = 0: This block has dimension 39, and it contains five sub-blocks of dimensions 
five, seven, eight (two of them), and eleven. We find eleven distinct eigenvalues. The 
simplest one is A ,i = q 4 — 9g 3 + 30g 2 — 44g + 25 with amplitude 2b^ = 2. Two 
eigenvalues are the solutions of the second-order equation 

C 2 -^(g 4 -10g 3 +39g 2 -68g+38)-g 7 +15g 6 -93g 5 +309g 4 -598g 3 +681g 2 -422g+109 = 

(3.11) 

with amplitude b^ = 1. The other eigenvalues come from polynomial equations of 
order three (with amplitude 2b^ = 2) and five (with amplitude = 1). 
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£ = 1: This block has dimension 494 and it contains eleven sub-blocks of dimensions 
ranging from 21 to 71. We find 48 different eigenvalues in this block. The simplest 
ones are given by A(i j2 ) = — (q 3 — Sq 2 + 20q — 14) and A(i i2 ,3) = — (q 3 — 6q 2 + 12q — 8), 
which also they appear in other blocks (see below). Two eigenvalues are the solutions 
of the second-order equation: 

^ 2 +^(2g 3 -16g 2 + 42g-38) + g 6 -16g 5 + 104g 4 -352g 3 + 656g 2 -642g+262 = (3.12) 

with amplitude = q — 1. The other eigenvalues come from polynomial equations of 
order ten, eleven (two of them), and twelve. 

£ = 1,2: The eigenvalue A(i i 2) = — (q 3 — 8q 2 + 20q — 14) appears in these two blocks 
with an amplitude (q - 2)(q + l)/2 = + bf\ 

£ = 2: This block has dimension 1432 and it contains 13 sub-blocks of dimensions 
between 56 and 180. We have found 120 distinct eigenvalues: A(i i 2), A(i 5 2,3), A(2,3) = 
g 2 — 5g+5 which also appears in the I = 3 sector (see below), and 117 eigenvalues coming 
from equations of order five (two equations), seven (two equations), nine, eleven, 13, 
14 (two equations), and 16 (two equations). 

£ = 1,2,3: The eigenvalue A(i 5 2,3) = — (g 3 — 6g 2 + 12g — 8) appears in these three 
blocks with amplitude (q - 3)(2g 2 - 3q + 4)/6 = + bf ] + b {3) . 

I = 2,3: The eigenvalue A( 2j 3) = q 2 — 5g + 5 appears in these two blocks with amplitude 
(2g 3 - 9g 2 + 7q - 6)/6 = bf ] + fo[ 3) . 

£ = 3: This block has dimension 644 and it contains seven sub-blocks of dimensions 
between 31 and 150. We have found 82 distinct eigenvalues coming from equations of 
order three, six, seven, eight (six equations), and nine (two equations). In addition, 
we obtain the already discussed eigenvalues A(i )2 ,3) and A( 2j 3). 

£ = 4: This block has dimension 163 and it contains four sub-blocks of dimensions 26, 
45, and 46 (two of them). We have found 49 distinct eigenvalues. The simplest one 
is A4 5 i = — (q — 2)(q — 3) with amplitude b^f 1 (2.33a). Six eigenvalues are solutions of 
three second-order equations 

= £ 2 _£(2g 2 -llg + 13) + g 4 -llg 3 + 43g 2 -71g + 40 (3.13a) 
= £ 2 _£( g 2 -7g + 8)-g 3 + 8g 2 -18g + ll (3.13b) 
= £ 2 _£( g 2 -7 g + io)-g 3 + 8g 2 -20g + 15 (3.13c) 

The other eigenvalues come from equations of order three (four equations), four (four 
equations), six, and eight. 

£ = 5: This block has dimension 16. The simplest eigenvalues are As^ = 5 — q (with 
amplitude bf^), X^^ = 4 — q (with amplitude 2b^), = 2 — q (with amplitude 
2&f' ) ), and A 5i4 — 1 — q (with amplitude bf^). Four eigenvalues are solutions of two 
second-order equations: 

= £ 2 + £(2g-7) + g 2 -7g + ll (3.14a) 
= £ 2 + £(2g - 5) + g 2 - 5g + 5 (3.14b) 
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The remaining eight eigenvalues come from two fourth-order equations. The amplitude 
of all these twelve eigenvalues is 2b^ (2.33e). 

• £ = 6: This one-dimensional block gives A 6 = 1 with amplitude (2.31e). 

For this strip we find 325 different eigenvalues; all of them but three (A(i )2 ), A(i ;2 ,3), and 
A(2,3)) belong to a single £ sector. 

In Figure 6(a) we have shown the chromatic zeros for the strips of lengths n = 30 and 
n = 60, and the limiting curve B 6 . This curve splits the complex g-plane into four regions. 
The region containing the real segment (0, 2) is dominated by an eigenvalue from the £ = 1 
sector. The region containing the real segment (2,5/2) is dominated by the £ = 2 sector. 
The two small complex conjugate regions are dominated by the £ = 1 sector. The rest of the 
complex g-plane is dominated by the £ = sector. 

The curve £> 5 crosses the real q axis at three points: q = 0, q = 2, and q « 2.9078023424. 
We also find four pairs of complex conjugate T points at q ~ 2. 1180081660 ± 2.2566257839 i, 
q « 2.1632612447±2.2209172289i,g « 2.8895175995±1. 1648906031 i, and q « 2.9197649745± 
0.7583704656 %. We also find a pair of complex-conjugate endpoints at q ~ 2.0571168133 ± 
2.3885607275 i. Finally, there is a single isolated limiting point at q — 1. 

3.6 L = 7 

The full transfer matrix T(7t) has dimension 20766. We have been unable to obtain the 
full eigenvalue structure, as some of the blocks for £ = 2,3 are very large. Thus, we cannot 
obtain the amplitudes. We find the following (partial) eigenvalue structure: 

• £ = 0: This block has dimension 111, and it contains six sub-blocks of dimensions 
six and 21 (five of them). We find 21 different eigenvalues coming from polynomial 
equations of order six and 15. 

• t = 1: This block has dimension 2794, and it contains 19 sub-blocks of dimensions 
103 and 196. We find 112 different eigenvalues coming from polynomial equations of 
order nine, 19, and 84. 

• £ = 2: This block has dimension 10224, and it contains 26 sub-blocks of dimensions 
270, 276, and 532. We find 292 different eigenvalues coming from polynomial equations 
of order 16, 22, 26, and 114 (two of them) arising from the characteristic polynomials 
of the sub-blocks of dimensions < 276. This description is not complete, as we cannot 
obtain the characteristic polynomial of the sub-blocks of dimension 532. However, 
we can numerically diagonalize this sub-block and find that it provides 12 additional 
distinct eigenvalues; hence, we expect to have 304 distinct eigenvalues in this sector. 

• £ = 3: This block has dimension 5615, and it contains 17 sub-blocks of dimensions 247 
and 483. We find 247 different eigenvalues coming from the polynomial equations of 
order 17, 23, 69, and 138 that arise from the characteristic polynomials of the smaller 
sub-blocks. Furthermore, we numerically find that the larger sub-blocks of dimension 
483 contains only six additional distinct eigenvalues. Thus, we expect 253 distinct 
eigenvalues for this sector. 
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• I = 4: This block has dimension 2711, and it contains eight sub-blocks of dimensions 
155, 158, and 308. We find 165 different eigenvalues coming from polynomial equations 
of order three, four, seven, eight, eleven, 33 (two of them), and 66. 

• I = 5: This block has dimension 288, and it contains four sub-blocks of dimension 
72. We find 72 different eigenvalues coming from polynomial equations of order four, 
eight, twelve, and 48. 

• I = 6: This block has dimension 22. We find four simple eigenvalues A6,i = q — 5, 
^6,2 = q — 4, A6,3 = q — 2, and A6,4 — q — 1. The other come from equations of order 
three and six (two of each). 

• t = 7: This one-dimensional block gives A 7 = —1. 

In Figure 6(b) we have shown the chromatic zeros for the strips of lengths n = 35 and 
n = 70, and the limiting curve B 7 . 

The curve Bj crosses the real q axis at three points: q = 0, q = 2, and q = 3. We 
also find six pairs of complex conjugate T points at q ~ 1.4876491726 ± 2.5577806355 i, 
q « 1.6498030004 ± 2.5002880707 i, q « 1.9765852319 ± 2.3664338035 i, q « 2.5531414480 ± 
1.8775702667 i, q « 2.8395564367 ± 1.4513813377^, and q « 2.8734078918 ± 1.3454749268 i. 
There are two complex conjugate oval-shaped regions between the T points at q w 1.488 ± 
2.558 « and q ~ 1.650 ± 2.500 i. And there is another pair of small bulb-like regions between 
the T points at q m 2.840 ± 1.451 % and q « 2.873 ± 1.345 i. 

Thus, the limiting curve B-j splits the complex g-plane into seven regions. The region 
containing the real segment (0, 2) is dominated by an eigenvalue from the £ = 1 sector. The 
region containing the real segment (2, 3) is dominated by the £ = 2 sector. The two small 
complex conjugate oval-shaped regions are dominated by the £ = sector; and the two small 
complex conjugate bulb-like regions are dominated by the £ = 1 sector. The rest of the 
complex g-plane is dominated by the £ = sector. 

We also find a pair of complex-conjugate endpoints at q ~ 1. 6947007027±2. 5327609879 i. 
Finally, there is a single isolated limiting point at q — 1. 

4 Triangular-lattice numerical results 

In this section we will analyze the results for the triangular-lattice strips of widths 2 < 
L < 7. We have checked our results using (in addition to the already known cases L < 4 
[16,18,19]) the identity (2.22). 

As for the square lattice, for each L we give a detailed description of the eigenvalues and 
amplitudes. For the already known cases L = 2, 3, we include for the sake of clarity the full 
expressions for the transfer matrices and the left and right vectors. The exact expressions 
for the new cases L = 5, 6 can be found in the mathematica file transf er6_triT.m available 
as part of the on-line version of this paper in the cond-mat archive at arXiv . org. As in 
the square-lattice case, for L = 7 we can only present a partial description of the eigenvalue 
structure, as some of the sub-blocks are very large. 
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4.1 L = 2 



The chromatic polynomial for this strip was computed by Chang and Shrock [19] (in fact, 
they computed the full Potts-model partition function). 

The full transfer matrix T(2t) has dimension five. In the basis {5i',2^2',i, ^i',i^2,2', ^i',2 + 
^2',i, #i',i +#2',2, 1}) w e find the following transfer matrix with a lower-triangular block struc- 
ture 
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q 2 — 5q + 6 / 



(4.1) 



The right and left u vectors are given by 



u T = (0,A,0,2A,A), with A = q(q - 1) 
= (0,1,0,0,0) 



(4.2a) 
(4.2b) 



The eigenvalues can be obtained from the diagonal blocks Ta. (each of them associated to 

n ot = i). 

After some algebra, we find the following eigenvalues and amplitudes: 

• £ = 0: The eigenvalue is Ao = q 2 — 5q + 6 with amplitude = 1. 

• £ = 1: There are two eigenvalues: Ai = — 2(q — 3) with amplitude = q — 1, and 
A Q = 0, which is common to the £ = 2 block. 

• £ = 1,2: The eigenvalue A Q = appears in these two blocks. Its amplitude is irrelevant 
for the computation of both the limiting curve £>2 and the chromatic polynomials 

P2pXn P (q)- 

• £ = 2: We find A = 0, and the eigenvalue A 2 = 2 with amplitude (2.32a). 

In Figure 7(a) we have shown the chromatic zeros for the strips of lengths n — 10 and 
n = 20, and the limiting curve Bi. This curve splits the complex g-plane into three regions. 
The outer one is dominated by the £ = sector. The one containing the real interval 
q E (0, 2) is dominated by the £ = 1 sector; and the one containing the real interval q G (2,4) 
is dominated by the £ = 2 sector. 

The curve B 2 crosses the real q axis at three points: q = 0, q = 2, and q = 4. This latter 
point is in fact a multiple point. Finally, there are two isolated limiting points at q — 1 and 
q = 3. 



4.2 L = 3 

The chromatic polynomial for this strip was first computed by Chang and Shrock [18]. 
The full transfer matrix T(3t) has dimension 13. In the basis {5i,3'52,i'^3,2' ? ^i,2'^2,3'^3,i', 

&L,l'$2,2'^3,3', ^1,2'^2,1' + • • • , 5l,2'^3,l' + • • • , £l,l'#2,2' + • • • , ^2,1'^3,2' + • • • , ^2,1'^3,3' + • • • , 5l,2'^3,3' + 
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. . . , 5i t 2' + . . . , <5i,i/ + . . • , #1,3' + . . . , 1} (where the dots ". . . " mean all possible connectivity 
states symmetric under rotations), it takes the form 
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where Sk is defined in (3.4a) and 



R = g 2 _7g+13 

A = q 3 - 9q 2 + 29q - 32 



(4.3) 

(4.4a) 
(4.4b) 



The right w id and left u vectors are given by 



u T = (0, 0, A, 0, 0, 3A, 0, 0, 0, 0, 3A, OA) , with A — q(q — l)(q - 2) (4.5a) 



w 



id 



(0,0, 1,0,0,0,0,0,0,0,0,0,0) 



(4.5b) 



The eigenvalues can be obtained from the diagonal blocks Tgj. (each of them associated to 
Vbot — 1)- We find the following eigenvalue structure: 

• £ = 0: There is a single eigenvalue Ao = q 3 — 9q 2 + 29q — 32 with amplitude = 1. 

• £ = 1: There are three distinct eigenvalues A^i = — 2(g — 4) 2 , and Ai ;2 ± = (q 2 — 5q + 
7) e ±2 ™/ 3 , with amplitude = g - 1. 

• £ = 2: The eigenvalues A( 2j3 ±) = e ±27 ™/ 3 also appear in the £ = 3 (see below). We 
also find A 2 ,i = g — 2, A 2j2 = 3g — 14, and A 2j3± = (7 — 2g)e ±27U//3 . The corresponding 
amplitudes are bf \ bf\ and b^\ respectively [c.f., (2.32)]. 

• £ = 2,3: The complex-conjugate pair A( 2j3± ) = e ±2m ^ appears in these two blocks 
with an amplitude q(q - l)(2q - 7)/6 = bf ] + ^ 3) 



^ + ■ Please note that in this 



case the linear combination for the amplitude is not unambiguously determined. 



• £ = 3: In addition to A( 2j3 ±), we find A3 = —2 with amplitude b\ (2.32c). 

In Figure 7(b) we have shown the chromatic zeros for the strips of lengths n — 15 and 
= 30, and the limiting curve B3. This curve splits the complex g-plane into three regions. 



n 



The region containing the real segment (0, 2) is dominated by an eigenvalue from the I 
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sector. The region containing the real segment (2, 3] is dominated by the i = 2 sector. The 
rest of the complex g-plane is dominated by the £ = eigenvalue A (4.4b). 

The curve B3 crosses the real q axis at three points: q = 0, q = 2, and q ~ 3.7240755514. 
We also find a pair of complex conjugate T points at q ~ 3.9777108575 ± 1.6284204463 i. 
Finally, there are two isolated limiting points at q — 1,3. 

4.3 L = 4 

The chromatic polynomial for this strip was computed by Chang and Shrock [16]. The 
full transfer matrix T(4 T ) has dimension 99. In this case, we restrict ourselves to report the 
eigenvalue structure: 

• £ > 0: The eigenvalue A Q = is common to all blocks < £ < 4. Its amplitude 
cannot be computed from our analysis; but it is unimportant to compute the chromatic 
polynomials and the limiting curve. 

• £ = 0: This block has dimension five with two sub-blocks of dimensions three and two. 
In addition to A Q = 0, we find two distinct eigenvalues coming from a second-order 
equation [16, Eq. (4.9)] and with amplitude = 1. 

• £ = 1: This block has dimension 26, and it contains three sub-blocks of dimensions six, 
and ten (two of them). We find eight different eigenvalues, one of them being A = 0. 
The other eigenvalues come from a third-order equation [16, Eq. (4.12)] and from a 
fourth-order equation [16, Eq. (4.14)]. The amplitudes are in all cases &W = q-l. 

• £ = 2: This block has dimension 52, and it contains three sub-blocks of dimensions 12, 
and 20 (two of them). There are 18 different eigenvalues. In addition to A Q = 0, we 
find two simple ones: A24 = 2(3 — q) and = 2(g — 3) 2 , both with amplitudes 
The other eigenvalues come from polynomial equations of order two, three, four, and 
six, respectively [16, Eqs. (4.11)/(4.13)/(4.15)/(4.17)]. The corresponding amplitudes 
are for the first three, and b^ for the last one [c.f. (2.32)]. 

• £ = 3: This block has dimension 12. We find three simple eigenvalues, A Q = 0, 
A.3,1 = 3 — q with amplitude 2&2 3 \ and A3 t 2 = 2(9 — 2q) with amplitude bf^ (2.32c). We 
also find six eigenvalues coming from polynomial equations of order two [16, Eq. (4.10)] 
and six [16, Eq. (4.16)], with common amplitude bf^ ■ The remaining two eigenvalues 
A(3,4±) are common to the £ = 4 block (see below). 

• £ = 3,4: There are two eigenvalues common to the £ = 3,4 sectors: A( 3)4 ±) = 1 ± % 
with amplitude (q - l)(q - 2)(3g 2 - llq - 6)/12 = b? + bf . 

• £ = 4: This is a four-dimensional block. In addition to A Q = and A (3^), we find 
A4 = 2 with amplitude b^ (2.33a). 

Thus, for this strip we find 38 distinct eigenvalues. All of them but three (A© and A(3 j4± )) 
belong to a single £ sector. 

In Figure 7(c) we have shown the chromatic zeros for the strips of lengths n = 20 and 
n = 40, and the limiting curve £>4. This curve splits the complex g-plane into three regions. 
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The region containing the real segment (0, 2) is dominated by an eigenvalue from the £ = 1 
sector. The region containing the real segment (2, 4) is dominated by the £ = 2 sector. The 
rest of the complex g-plane is dominated by the £ = sector. 

The curve £> 4 crosses the real q axis at three points: q = 0, q = 2, and q = 4. We also 
find a pair of complex conjugate T points at q « 3.3943280448 ± 2. 1041744504 %. Finally, 
there are two isolated limiting points at q — 1, 3. 

4.4 L = 5 

The full transfer matrix T(5 T ) has dimension 653. The relevant eigenvalues can be 
extracted from the diagonal blocks T U (5 T ). We find the following eigenvalue structure: 

• £ = 0: This block has dimension eight, and it contains two sub-blocks of dimensions 
two and six. We find six distinct eigenvalues (with amplitude = 1) coming from a 
fourth-order equation and the following second-order equation: 

= £ 2 _£( g 5 _i5g 4 + 98g 3 - 355g 2 + 711g-614) + 3g 8 -75g 7 + 817g 6 

-5062g 5 + 19492g 4 - 47696g 3 + 72273g 2 - 61822g + 22764 (4.6) 

• t = 1: This block has dimension 125 and it contains seven sub-blocks of dimension 25. 
We find 24 different eigenvalues in this block. The eigenvalue \n,2) — —{q 2 — 5q+6) also 
appears in the £ = 2 block (see below). The other eigenvalues come from polynomial 
equations of order three and 20. They all have the same amplitude U 1 ^ = q — 1. 

• £ = 1,2: The eigenvalue A(i j2 ) = — {q 2 — 5g + 6) with amplitude 2q(q — 2) = 2(&W + 
frj 2 ^ + b^) appears in these two blocks. 

• £ = 2: This block has dimension 360 and it contains six sub-blocks of dimension 60. 
We find 57 distinct eigenvalues, including \i,2)- We also find the simple eigenvalue 
A2,i = — (q — 4)(2g — 7) with amplitude (2.32b). The other eigenvalues come from 
polynomial equations of order three, four, and 24 (two of them) . 

• £ = 3: This block has dimension 135 and it contains three sub-blocks of dimension 45. 
The 42 distinct eigenvalues come from the solutions of polynomial equations of order: 
three (two of them), 12, and 24. 

• £ = 4: This block has dimension 20. Three eigenvalues are rather simple: A41 = q — 4 
with amplitude Ib^' ', = q — 2 with amplitude b^\ and A^ = 5q— 22 with amplitude 
b^ . Four eigenvalues A^^j) are common to the £ = 5 block (see below). The other 
eigenvalues come from polynomial equations of order four and eight, with amplitudes 
b^ and b^\ respectively [c.f., (2.33)]. 

• £ = 4, 5: The four eigenvalues A(4 ;5) j) are the solution of the equation 

£ 4 + 3£ 3 + 4£ 2 + 2£ + 1 = (4.7) 

Their amplitude is q(q - 3)(g - 2)(4g 2 - 15q + l)/20 = b^ + bf = bf + bf\ Again, 
we find an ambiguity in the linear combination of basic amplitudes. 
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• £ = 5: This block is five-dimensional. In addition to A(4 i5; j), we find A5 = —2 with 
amplitude bf (2.33d). 

Thus, for this strip we find 148 different eigenvalues. All of them but five (A(i ; 2) = 
— (q 2 — 5q + 6) and A( 4i5i j)) belong to a single £ sector. 

In Figure 7(d) we have shown the chromatic zeros for the strips of lengths n = 25 and 
n = 50, and the limiting curve B$. This curve splits the complex g-plane into six regions. 
The region containing the real segment (0, 2) is dominated by an eigenvalue from the £ = 1 
sector. The regions containing the real segment (2, 4) is dominated by the £ = 2 sector. The 
rest of the complex g-plane (including the two complex-conjugate triangular-shaped regions) 
is dominated by the £ = sector. 

The curve B 5 crosses the real q axis at four points: q = 0, q = 2, q 3.8482632901 
and q — 4. We also find four pairs of complex conjugate T points at q ~ 2.7382976052 ± 
2.8142366786 i, q « 2.7643522733 ± 2.3252257137i, q « 3.7492171129 ± 1.9879721106 i, and 
q ~ 3.9634199872±0.3751435740i. Finally, there are two isolated limiting points at q — 1, 3. 

4.5 L = 6 

The full transfer matrix T(6t) has dimension 5194. The relevant eigenvalues can be 
extracted from the diagonal blocks T«(6t). 

• £ > 0: There is a single eigenvalue common to all blocks A Q = 0. As for the cases 
L — 2, 4, its amplitude cannot be determined; but it is irrelevant to compute both the 
chromatic polynomials and the limiting curve. 

• £ = 0: This block has dimension 49, and it contains five sub-blocks of dimensions 
ranging from five to 15. We find 14 distinct eigenvalues: A = 0, the solutions of the 
second-order equation 

= £ 2 - £(3g 4 - 33g 3 + 132g 2 - 225g + 135) + 3g 8 - 66g 7 + 627g 6 

-3354g 5 + 11028g 4 - 22770g 3 + 28755g 2 - 20250g + 6075 (4.8) 

and the solutions of equations of order five and six. All these eigenvalues have the 
amplitude 6^ = 1. 

• £ = 1: This block has dimension 843 and it contains 13 sub-blocks of dimensions 
ranging from 25 to 71. We find 62 different eigenvalues in this block: A = 0, the 
solutions of the second-order equation 

^ 2 + ^(3g 3 -27g 2 + 79g-75)+g 6 -18g 5 + 131g 4 -492g 3 + 1000g 2 -1032g + 414 = (4.9) 

and the solutions of equations of order four, nine, 22, and 24. 

• £ — 2: This block has dimension 2712 and it contains 21 sub-blocks of dimensions 
ranging from 96 to 180. We find 165 different eigenvalues in this block: A = 0, and 
the solutions of polynomial equations of order three (two of them), four (two of them), 
ten (three of them), 28 (two of them), and 32 (two of them). 
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• i = 3: This block has dimension 1254 and it contains nine sub-blocks of dimensions 
ranging from 54 to 150. We find 136 different eigenvalues in this block: A Q = 0, 
^3,1 = (q — 2)((? — 3) with amplitude 56^ , and solutions of equations of order four, six, 
eight, 16 (five of them), and 18 (two of them). 

• t = 4: This block has dimension 300 and it contains five sub-blocks of dimensions 
ranging from 48 to 84. We find 76 different eigenvalues in this block: A = 0, A^i = 
3q 2 — 18q + 26 with amplitude (2.33c), A^ = 2(3 — q) with amplitude 3b^\ the 
solutions of the second-order equation 

£ 2 + 4g 4 - 44g 3 + 176g 2 - 304g + 192 = (4.10) 

with amplitude and the solutions of equations of order three, four, six (two of 
them), eight (three of them), twelve, and 16. 

• £ = 5: This block has dimension 30, and it contains 28 different eigenvalues: A = 0, 
A 5j i = 26 — 6q with amplitude bf^ (2.33d), the solutions of the second-order equation 

£ 2 + £(2g-7) + g 2 -7g + ll = (4.11) 

with amplitude 2b^\ four eigenvalues A^,?) common to the £ = 6 sector (see below), 
and the solutions of equations of order four and eight (two of them). 

• £ = 5,6: There are four common non-zero eigenvalues to these two sectors, coming 
from two second-order equations: A(5 ; 6,i±) = V / 3e ±Mr / 6 with amplitude (q — 3)(5g 5 — 
39g 4 + 93g 3 - 76g 2 + 22g + 10), and A( 5)6i 2±) = e ±iw / 3 with amplitude (q - l)(5g 5 - 
49g 4 + 161g 3 - 184g 2 + 6q + 30). 

• t = 6: This block has dimension 6. The eigenvalues are A Q = 0, A^^j), and \$ = 2 
with amplitude q(q — 3)(q — l)(g 3 — 8g + 19g — 1 1)/6. 

Thus, in this strip we find 476 distinct eigenvalues. All of them but five (A = and 
A(5,6,i)) belong to a single sector. 

In Figure 8(a) we have shown the chromatic zeros for the strips of lengths n = 30 
and n = 60, and the limiting curve £> 6 . This curve splits the complex g-plane into seven 
regions. The region containing the real segment (0, 2) is dominated by an eigenvalue from 
the t = 1 sector. The region containing the real segment (2, 3) is dominated by the £ = 2 
sector. The two bulb-like complex-conjugate regions appearing on the right hand-side of the 
limiting curve belong to the £ = 1 sector and they are rather special (see below). The rest 
of the complex g-plane (including the two complex-conjugate triangular-shaped regions) is 
dominated by the £ = sector. 

The curve Bq crosses the real q axis at three points: q = 0, q = 2, and q ~ 3.7943886378. 
We also find five pairs of complex conjugate T points at q ~ 2.0352758561 ± 3.1492288108 i, 
q w 3.3399899103 ± 2.4428522934 i, q « 2.5037195029 ± 2.3940503612 i, q « 4.1131475777 ± 
0.7636584371 i, and q w 4.1295836790 ± 0.8010545279 i. Finally, there are two isolated 
limiting points at q — 1, 3. 

Important technical remark. As stated above, the two bulb-like complex-conjugate 
protruding from T points q ss 4.11 ± 0.76 z and q ps 4.13 ± 0.80 z, are rather special. Their 



30 



peculiarity arises from the fact that inside these two bulb-like regions there is a pair of 
exactly equimodular dominant eigenvalues [with non-zero amplitudes ~ (q — 1)]. Indeed, we 
do not see chromatic zeros becoming dense in the interior of these regions, but only on their 
boundaries. 

The explanation of this phenomenon is rather simple: the two dominant eigenvalues 
come from a polynomial equation of order 24 belonging to the £ = 1 sector. A closer look 
at this equation reveals that its roots Zk can be written clS Zfc — Zq /j. 

e ±vn/3 w h ere the z k are 

the roots of another polynomial equation of order 12 with integer coefficients (that can be 
easily computed). It is important to stress that the existence of this pair of equimodular 
eigenvalues does not mean that the chromatic zeros accumulate in this open set in the 
infinite-length limit. According to the Beraha-Kahane- Weiss theorem [25], the eigenvalues 
should satisfy a "no-degenerate-dominance" condition: there do not exist indexes k ^ k' 
such that Afc = ooXk' with \u\ = 1, and such that the region where A& and A&/ are dominant 
has non-empty interior. This is precisely the condition violated by the leading eigenvalues 
in the interior of the above bulb-like regions. Thus, in order to compute the limiting curve 
in these regions, we should only consider the roots Zo^k of the 12-order equation. In this way, 
we ensure the "no-degenerate-dominance" condition. In conclusion, the chromatic zeros 
cannot accumulate in the interior of these bulk-like regions; but on their boundary (where 
the dominant z 0t k r °ot from the £ = 1 sector becomes equimodular to the leading eigenvalue 
from the £ = sector) . 

In other strips we have already found the existence of two different eigenvalues differing 
by a phase [i.e., A and Xe 1 ^}. However, this is the first case where one of these "degenerate" 
pairs is dominant on an open set of the complex g-plane. In order to correctly compute the 
limiting curve, one should eliminate one of eigenvalues of each "degenerate" pair to ensure 
the hypothesis of the Beraha-Kahane- Weiss theorem. 

4.6 L = 7 

The full transfer matrix T(7t) has dimension 40996. 10 We have been unable to obtain 
the full eigenvalue structure, as some of the blocks for £ = 2,3 are very large. We find the 
following (partial) eigenvalue structure: 

• t = 0: This block has dimension 186, and contains five sub-blocks of dimension 36, 
and a single block of dimension six. We find 36 distinct eigenvalues coming from two 
polynomial equations of order six and 30, respectively. 

• t = 1: This block has dimension 5488, and contains 28 sub-blocks of dimension 196. 
We find 196 distinct eigenvalues coming from a single polynomial equation of order 
196. 

• t = 2: This block has dimension 20216, and contains 38 sub-blocks of dimension 532. 
A numerical study of these sub-blocks reveals that there are 508 distinct eigenvalues 
in this sector. 



10 This matrix is so large that its mathematica representation cannot be fitted in 12Gb of RAM. In order 
to deal with it, we have used mathematica's SparseMatrix function, which saves a lot of memory due to 
the block-diagonal structure of our matrix. 
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• i = 3: The dimension of this block is 11109 and it contains 23 sub-blocks of dimensions 
483. We numerically find that there are 454 distinct eigenvalues for this sector. 

• i = 4: This block has dimension 3388, and contains eleven sub-blocks of dimension 308. 
There are 289 distinct eigenvalues: two of them are rather simple: A 4j i = — (q 2 — 5g + 6) 
and A 4 i = — (6g 2 — 44g + 80), while the rest come from polynomial equations of order 
three, four, five, eleven, 66 (two of them), and 132. 

• I = 5: The dimension of this block is 560 and it contains four sub-blocks of dimensions 
140. There are 122 distinct eigenvalues coming from polynomial equations of order four, 
eight, 24, and 96. 

• £ = 6: This block has dimension 42. We find four simple eigenvalues A6,i = q — 5, 
^6,2 — Q — 3, \s t 3 = q — 2, and Xq^ = 7q — 30. There are 30 additional eigenvalues 
coming from equations of order six and twelve (two of them). Finally, there are six 
other eigenvalues common to the £ = 7 block. 

• I = 6,7: We find six eigenvalues common to these blocks. They are the solutions of 
the equation 

£ 6 + 5£ 5 + ll£ 4 + 13£ 3 + 9£ 2 + 3£ + l = (4.12) 

• £ = 7: This block has dimension seven, and in addition to the eigenvalues common to 
the £ = 6 block, we find A 7 = — 2. 

In Figure 8(b) we have shown the chromatic zeros for the strips of lengths n = 35 and 
n = 70, and the limiting curve B-j. This curve splits the complex g-plane into eight regions. 
The region containing the real segment (0, 2) is dominated by an eigenvalue from the i = 1 
sector; the region containing the real segment (2, 3) is dominated by the £ = 2 sector; and the 
small region containing the real segment (3.8, 4) is dominated also by the £ = 2 sector. There 
are two small lens-like complex-conjugate regions which belong to the £ = 1 sector. The rest 
of the complex g-plane (including the two complex-conjugate triangular-shaped regions) is 
dominated by the £ = sector. 

The curve £>7 crosses the real q axis at four points: q = 0, q = 2, q 3.8605508375, 
and q — 4. We also find six pairs of complex conjugate T points at q 1.4835869410 ± 
3.2225993236 i, q w 2.4516540074 ± 2.4976908148 i, q w 3.0454262812 ± 2.6881696318 i, q w 
3.3660354065 ± 2.4484944601 i, q w 3.8749646144 ± 1.8155509876 i, and q w 3.9442064557 ± 
0.2720862591 2. Finally, there are two isolated limiting points at q — 1,3. 

5 Discussion 

5.1 Square lattice limiting curves 

A summary of qualitative results for the limiting curves £>, as obtained in Section 3, is 
shown in Table 6. In particular, we have computed qo{L), the maximum real value in B for 
strip widths L < 7 from the symbolic transfer matrices. In Figure 9, we display together 
all the computed limiting curves Bl with 2 < L < 7. For comparison, we also include the 
limiting curve for a strip of width L = 10 with cylindrical boundary conditions [12]. 
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These computations can be taken a lot further by diagonalizing the diagonal blocks 
T^(Lt) of the numerical transfer matrix. Indeed it suffices to diagonalize one of its sub- 
blocks, e.g., the one with a trivial bottom- row connectivity. This amounts to working with a 
transfer matrix that acts on connectivities of only L (rather than 2L) points, with exactly I 
marked clusters. For the purpose of getting only the eigenvalues, without going outside the 
specified sub-block, transfer matrix elements which would amount to adjoining two distinct 
marked clusters are set to zero. 

Since further we need only the largest (or sometimes the first few largest) eigenvalues 
of the given sub-block, we can diagonalize using a standard iterative procedure (the so- 
called power method [40]) which has the immense advantage of allowing for sparse matrix 
factorization (i.e., the lattice edges are added one by one) combined with standard hashing 
techniques. This allowed us to access widths L < 12. In all cases, qo(L) was found to 
correspond to the crossing (in modulus) of the leading eigenvalues of the I = and 1 = 1 
sectors. The precise value of qo(L) was then narrowed in by a Newton- Ralphson method. 
The final results are shown in Table 6. 11 

From the curious fact that qo(L) = 3 exactly for L — 3, 5, 7, 9, 11 we conjecture that this 
is true for all square-lattice strips with odd widths: 

Conjecture 5.1 For square-lattice strips with fully periodic boundary conditions of widths 
L = 2p + 1 with p > 1, we have qo(L) = 3. 

In other words, we have found a particular sequence of strip graphs for which the value of 
qo(L) is constant and, furthermore, equal to the expected critical value <? c (sq) = 3. Indeed, 
recall that for q = 3, the square-lattice chromatic polynomial is equal to three times the 
partition function of the six- vertex model with all vertex weights equal to one [33]; it is 
well-known that this latter model is critical (with central charge c = 1). 

We next conjecture that parity effects do not affect the value of q in the thermodynamic 
limit: 

Conjecture 5.2 For square-lattice strips with fully periodic boundary conditions of widths 
L = 2p, the limit go(sq) = lim qo(L) exists and is equal to g c (sq) = 3. 

This conjecture is very similar to the one made in [12]; the only difference being the 
boundary conditions, cylindrical rather than toroidal. 

As to the value of qo(L) for even L, the data of Table 6 is well fitted by the power-law 
form q (L) = q c (sq) + AL~ l / u . Our estimates are g„(sq) = 2.999 ±0.012 and v = 0.49 ±0.02. 
These results strongly corroborate Conjecture 5.2. A better estimate for v can be obtained 
by fixing go(sq) = 3 in the above Ansatz: we obtain the critical exponent v = 0.502 ± 0.002. 
We conjecture that the exact value of the exponent is v — 1/2, consistent with the transition 
being first-order in q. 

The limiting curves El for 2 < L < 7 show some regularities (see Figure 9): for real q, we 
have three different phases and each of them corresponds to a different £ sector (£ = 0, 1, 2). 
Furthermore, within each phase, the leading amplitude does not depend on L. All our 
empirical findings can be summarized in the following conjecture: 

11 Indeed, for 2 < L < 7 the numerical values found using the numerical transfer matrices agree with those 
found using the symbolic matrices. 
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Conjecture 5.3 The phase diagram of the zero-temperature square-lattice Potts antiferro- 
magnet with toroidal boundary conditions at real q has three different phases characterized 
by the number of bridges £ as follows: 

(a) q E (—00, 0) U (q (L), 00) characterized by £ = with amplitude = 1 

(b) q E (0, 2) characterized by £ = 1 with amplitude = q — 1 

(c) q E (2, qo(L)) characterized by £ = 2 with amplitude bf^ = q(q — 3)/2 

Note that for all widths studied, 2 < L < 12, we have that 2 < qo(L) < 3. 

We also find a curious behavior of one of the eigenvalues belonging to the highest sector 
for each strip width: 

Conjecture 5.4 The diagonal sub-block Tl^l for a square-lattice strip with fully toroidal 
boundary conditions and width L contains the eigenvalue A = (— 1) L with amplitude b^ L \ 

5.2 Triangular lattice limiting curves 

A summary of qualitative results for the limiting curves B is shown in Table 7. In 
Figure 10, we display together all the computed limiting curves Bl with 2 < L < 7. For 
comparison, we also include the infinite- volume limiting curve obtained by Baxter [14]. The 
results in Table 7 with L < 7 were obtained from the symbolic transfer matrices in Section 4, 
and those with L < 12 were found by the numerical procedure outlined in Section 5.1. In 
all cases qo(L) is given by the crossing between the dominant eigenvalues in the £ = and 
£ = 2 sectors. 

From the curious fact that qo(L) = 4 exactly for L = 2,4,5,7,8,10,11 we conjecture 
that: 

Conjecture 5.5 For triangular-lattice strips with fully periodic boundary conditions, and 
for all widths L = 3p ± 1 with p > 1, we have qo(L) = 4. 

A quick glance at Table 7 however reveals that the values of qo(L = 3p) are definitely 
not converging to 4. We therefore examine the data with L = 3p more closely. A first fit 
with q (L) = q (tn) + AL,- l l v yields g (tri) ~ 3.820 and v ~ 0.52 (See Table 8). This is 
consistent with the expected value v — 1/2 corresponding to a first-order phase transition 
(in the variable q), as we have already found for the square lattice. Fixing 1/ = 1/2 we then 
arrive at the final estimate go(tri) = 3.8196 ± 0.0005. 

Remark. We have performed all the fits reported in this paper using the following method: 
as the data is essentially exact, we have done a standard least-squares fit with as many 
data points as the number of unknown parameters in the Ansatz, so that in each fit 
there are no effective degrees of freedom. To estimate the error bars (due to the existence 
of correction-to-scaling terms not included in the Ansatz), we introduce the parameter L min 
such that each fit is performed with Na consecutive data points with L > L min . From the 
variation of the estimates as a function of L min , we deduce the corresponding error bars. 
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Finally, from the FSS theory of first-order phase transitions [24], we expect that higher 
corrections-to-scaling terms should be integer powers of L~ l l u = L~ 2 . Thus, we have fitted 
our data to the improved Ansatz 

q (L) = g (tri) + AL~ 2 + BL~ A (5.1) 

The results are displayed in Table 8, and from these we conclude that g (tri) = 3.8198 ± 
0.0006. 

We now recall that the triangular-lattice chromatic polynomial was actually shown by 
Baxter [14] to be an integrable model. In particular, using coordinate Bethe Ansatz, this 
author identified three different analytic expressions for the free energy in the thermodynamic 
limit. 

To be precise, Baxter used first the Nienhuis mapping [41] between the zero-temperature 
g-state Potts antiferromagnet on the triangular lattice and the critical O(n) loop model on 
the hexagonal lattice. That mapping includes an exact relation between q and the weight n 
of a (bulk) loop. Baxter then transformed the loop model into a vertex model, by the usual 
trick of redistributing the loop weight n as local complex vertex weights. Unfortunately, 
when imposing periodic transverse boundary conditions on the vertex model, the loops of 
non-trivial homotopy (i.e., that wrap around the periodic direction) get the incorrect weight 
of 2, a fact not mentioned in Baxter's paper. Actually, while giving a (boundary) weight of 
n' — n to those loops (by twisting the vertex model) would have been a nicer way to perform 
the analysis of the O(n) model, the subtleties of the Nienhuis mapping [41] are such that yet 
a different weight n'(q) ^ n would have been needed to ensure that clusters of non-trivial 
homotopy in the Potts antiferromagnet carry their correct weight of q. 

One might think that these subtleties have no implications on the behavior in the ther- 
modynamic limit, but the reader is reminded that we are here dealing with the antiferromag- 
netic regime in which boundary conditions can (and sometimes do) matter. Nevertheless, it 
seems reasonable to expect that this should have no consequence for the bulk free energies. 
Numerical analysis [13, 14] of Baxter's analytically determined free energies shows that the 
corresponding eigenvalues become equimodular at 

q B ~ 3.8196717312. (5.2) 

The value (5.2) is determined implicitly by equating two analytical expressions given in 
Ref. [14]. Given the excellent agreement with our above numerical determination of g (tri) 
we can therefore conjecture that: 

Conjecture 5.6 For triangular-lattice strips with fully periodic boundary conditions of widths 
L = 3p, the limit g (tri) = lim qo(L) exists and is equal to qs of Eq. (5.2). Furthermore, 

p— »oo 

the sequence qo(3p) is monotonically increasing in p. 

The number of regions enclosed by B having a non-zero intersection with the real g-axis 
is 2 for L = 2,3,4 [cf. Figure 7 and Table 7]. However, for L = 5 a further such region 
appears, whose intersecting with the real g-axis is the interval (52 (L), qo(L)), where q-i (L) ~ 
3.84826329. The value q2{L) corresponds to the crossing of the two largest eigenvalues in the 
£ = 2 sector and is tabulated for strip widths L < 12 in Table 7. This table gives compelling 
evidence that q2(L) tends to a limit, ^(tri) < 4, independently of the parity of Lmod3. 
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As above, we first fit the data for q2(L) with a fixed parity, Lmod3 = 0, to the form 
q2(L) = g 2 (tri) + AL^ 1 ^ (See Table 8). This is again favourable of v ?s 1/2. Repeating 
the fit with a fixed v — 1/2 we arrive at the final estimate (fe(tri) = 3.8194 ± 0.0016. Using 
the improved Ansatz (5.1), we get the value (^(tri) = 3.8199 ± 0.0017. We can therefore 
conjecture that: 

Conjecture 5.7 For triangular-lattice strips with fully periodic boundary conditions the 
limit (^(tri) = lim q 2 {L) exists and is equal to qs of Eq. (5.2). Furthermore, the sequence 

<?2 (3p) is monotonically decreasing in p. 

The monotonic nature of the sequences go(tri) and q 2 (tn), and the existence of a common 
limit (as stated in Conjectures 5.6-5.7), means that this common limit is bounded by the 
finite-widths values: i.e., qo(3p) < q^ < q2(3p) for all p > 2. The data for L = 12 displayed 
in Table 7 implies that 

3.81314 < q B < 3.82428 (5.3) 

In particular, this narrow interval excludes that this common limit would be equal to the 
value qc = Bi 2 (1-12) proposed in Ref. [9]. 

As for the square-lattice case, the limiting curves B L for 2 < L < 7 show some regularities 
(see Figure 10): for real q, we have up to four different phases belonging to only three different 
£ sectors [I = 0, 1,2). Again, within each phase, the leading amplitude does not depend on 
L. All our empirical findings can be summarized in the following conjecture: 

Conjecture 5.8 The phase diagram of the zero-temperature triangular-lattice Potts antifer- 
romagnet with toroidal boundary conditions at real q has four different phases characterized 
by the number of bridges £ as follows: 

(a) q E (— oo, 0) U (q (L), oo) characterized by i = with amplitude = 1 

(b) q E (0, 2) characterized by £ = 1 with amplitude b^ = q — 1 

(c) q E (2, q2{L)) characterized by I = 2 with amplitude = q(q — 3)/2 

(d) q E (q 2 (L), qo(L)) characterized by £ = 2 with amplitude bf^ = q(q — 3)/2 

Note that the phase (d) only appears for finite width L whenever qo(L) > q2{L) (e.g. L = 
5, 7). But, according to Conjectures 5.5-5.7, this phase does survive in the thermodynamic 
limit L — > oo. We believe that this phase will appear for all L > 5 such that Lmod3 ^ 0. 

For 5 < L < 7 we further note the appearance of two disjoint triangular-shaped complex- 
conjugate enclosed regions, whose boundary contains three T points, one of which belongs 
to the branch of B that contains the point q = 2. The dominant eigenvalue inside these 
enclosed regions comes from the £ = sector. Figures 7-8 provide some evidence that these 
regions may subsist in the thermodynamic limit. 

Finally, for L = 6 we note the emergence of a pair of complex-conjugate enclosed regions 
analogous to those found for the square lattice. More precisely, these enclosed regions contain 
points for which Re q > 4, and the corresponding dominant eigenvalue comes from the £ = 1 
sector. We find it likely that such regions will appear for all L = 3p with p > 2, and 
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will merge in the L — > oo limit so as to form a branch that intersects the real g-axis in 
q = q (tn) = 4. 

There is a analogous conjecture for the triangular lattice to Conjecture 5.4 about one of 
the eigenvalues of the highest sector: 

Conjecture 5.9 The diagonal sub-block T L L for a triangular-lattice strip with fully toroidal 
boundary conditions and width L contains the eigenvalue X = (— 1) L 2 with amplitude b[ L \ 



5.3 Triangular-lattice free energy 

In this section we will discuss the free energy for the triangular-lattice strips we have 
considered and its relation to the analytic eigenvalues found by Baxter [14]. 

The (real part of the) free energy per unit area fii^q) for a triangular-lattice strip of width 
L is given by 12 

h(q) = log|A*(L;(z)| (5.4) 

where A*(L; q) is the dominant eigenvalue of the transfer matrix T(L), and G is the geometric 
factor [c.f., [9,10]] 

G(tri) = ^ (5.5) 

relating the free energy per unit area and per spin. With this normalization the FSS of the 
free energy is given by 

c(q)n 
6L 2 

where /bulk is the bulk free energy per unit area (which is expected to be independent of 
boundary conditions), c(q) is the (effective) central charge, and the dots stand for higher- 
order corrections depending on the value of q (in particular, CFT predicts a non-universal 
L~ 4 correction term). 

In Figure 11, we have plotted the two most dominant free energies for each of the three 
first sectors £ = 0, 1, 2 for a triangular-lattice strip of width L = 12 around q = qB (1.11). 13 
The other sectors £ > 3 are always sub-dominant and play no role in this discussion. The 
relevant eigenvalue crossing occurs at q = q (L = 12) « 3.8131463701 < ge where the 
dominant eigenvalue changes from the £ = 2 sector (solid black line) to the £ = sector (dot- 
dashed red line). There is also a crossing between the two most dominant £ = 2 eigenvalues 
at q — q2{L = 12) « 3.82428323 > qB- This crossing gives rise to a new £ = 2 phase (between 
q and g 2 ) for all strips of width not a multiple of three. In both cases, the £ = 1 sector plays 
no role in the phase structure around q = qB- This computation gives numerical evidence 
that our scenario is correct and no new elements are at play. 



h(q) = /buik(g) + ^tt + ■ ■ ■ ( 5 - 6 ) 



12 Note that in previous papers [9, 11-13], we have considered the free energy per spin and used the same 
notation /l(<?)- 

13 The width L = 12 is the largest one for which this computation could be carried out numerically. To 
keep computation times at a reasonable level wc therefore used a rather coarse mesh, Aq = 0.01, for scanning 
through the values of q. This is the reason why the curves in Figure 11 docs not look as smooth as they 
should do around the eigenvalue-crossings within a given level I. This is particularly clear for 1 = 0. 
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Another interesting question is whether the bulk free energy /buik(<?) agrees or not with 
Baxter's free energies per unit area (1/G) log |</j| (i = 1, 2, 3) [14]. (See Ref. [13, Section 6.1] 
for a detailed description of Baxter's solution). In particular, Baxter found that the free 
energy is given by in the region containing the interval q G (0,gB), by g<i in the region 
containing the interval q G (ge,4), and by g\ in the rest of the complex g-plane. 

In Figure 11 we observe that the agreement with Baxter's solution is very good: for q > q-Q 
the free energy is very similar to the g<i contribution (actually it seems to be almost equal to 
the contribution of the £ = 1 sector), while for q < q^, it is similar to the contribution of g%. 

A striking difference between our phase-transition diagrams and Baxter's solution is the 
transition line going through q = 2, which is absent in the latter. This line separates the 
phase characterized by £ = 1 from the one characterized by £ = 2 (which agrees with Baxter's 
solution close to q — ge)- Therefore a natural question is whether there is one eigenvalue 
missing in Baxter's solution and corresponding to the £ = 1 sector. To answer this question 
we have plotted in Figure 12(a) the ratio of the free energy per unit area /l(<?) and Baxter's 
free energy (1/G) log (73(g) in the interval q G (0,3.5). To avoid parity effects, we have 
considered only lattices of widths L = 3p with p = 2,3,4,5, and for each width we have 
considered the free-energy contribution from the three lowest sectors £ = 0, 1, 2 (represented 
by solid, dashed, and dot-dashed curves, respectively). 

Note that the eigenvalue crossings occur at the points q — 0,2. Another exact crossing, 
involving the second and third most dominant eigenvalues, takes place at q — 1. The fact 
that the loci of these crossings are subject to no finite-size corrections can presumably be 
explained using quantum group arguments. 

The numerical agreement with g 3 is not very good as the ratio is consistently below the 
expected result (i.e., 1). The reason is that the FSS corrections (5.6) have not been taken 
into account. A better agreement can be obtained if we take into account the first correction 
term in (5.6). We then define the corrected free energy per unit area 

/«t,l(?) = h(q) - ^ff^ = /bui k (g) + o (L- 2 ) (5.7) 

where c e g(q) is the effective central charge for the BK phase discussed in the next section; see 
(5.14). We expect that this value will be closer to the thermodynamic limit /bulk (9)5 as the 
correction terms are o(L~ 2 ). This plot is shown in Figure 12(b). Our first observation is that 
the values of the ratio are much closer to the expected value of 1 than in Figure 12(a). This 
shows that our assumption that CFT holds for this model is valid. Note that the eigenvalue 
crossings do not need to be at the expected points q — 0, 2, as the correction in (5.7) depends 
on both q and L. The second observation is that the agreement of the £ = 2 free energy with 
the contribution of g 3 is very good for q > 1 (less than 0.02% for L — 9, 12). For q ^ 1, the 
best agreement corresponds to the £ = sector. However, it is difficult from this data alone 
to draw any definitive conclusion about the necessity of including an extra eigenvalue to the 
set introduced by Baxter. 

5.4 Conformal properties 

Following Refs. [29,30] we expect a number of the conformal properties of the chromatic 
polynomial to be given by the BK antiferromagnetic phase. In brief, the situation is the 
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following. 14 The two-dimensional Potts model is exactly solvable (i.e., integrable) on the 
curves 

v = ± A /o (square lattice) 

V V ' (5.8) 

v 3 + 3v 2 = q (triangular lattice) 

The part of the curves with < q < 4 and v > corresponds to the usual second- 
order transition between a ferromagnetic and a paramagnetic phase. Its critical properties 
are given by the Coulomb gas (free field) construction with coupling constant g G [0, 1/2]; 
in particular the thermal operator (corresponding to a perturbation of the critical point in 
the v variable) is relevant in the RG sense. Now, the analytic continuation to the region 
with g G (1/2, 1) describes the lower (resp. middle) branch of the curves (5.8) for the square 
(resp. triangular) lattice. Crucially, at these critical points the thermal operator is irrelevant. 
Therefore, for each value of q G (0,4), these critical points will act as an RG attractor for 
a finite interval (v -(g), v + (q)): this is the BK phase. For the square lattice, the end-points 
of these intervals are given by the two mutually dual antiferromagnetic phase transition 
points [31] 

v ± (q) = -2±^~q (5.9) 

whereas for the triangular lattice v±(q) are presently unknown. 

The upshot is that the part of the chromatic line v — — 1 that intersects the BK phase 
will have its critical properties determined by the latter. For the square lattice this happens 
for Reg G (0,3), according to Eq. (5.9), in perfect agreement with Conjecture 5.2. For the 
triangular lattice, we obtain a consistent scenario if we suppose that the intersection with 
the BK phase is for Reg G (0, g2), with g 2 given by Conjecture 5.7. 

Within the BK phase, we can parametrize q by 

.2 



q = 4 cos 2 (5.10) 
where p G (2, oo). The central charge is then 

c(p) = 1 - e ±^l. (5.H) 

V 

Furthermore, the theory possesses an infinite set of bulk operators Ok, with k — 1, 2, . . ., of 
conformal weight [42] 

hM = k2 ' { l~ l? - (5-12) 

Geometrically, the insertion of an Ok operator at either end of the strip corresponds to a 
situation in which 2k cluster boundaries propagate along the length direction of the strip. 
(The boundaries may of course also wind around the periodic transverse direction.) For 
k > 2 we may equivalently think of k propagating clusters. It should be carefully noted 
that the case k — 1 is special: a single propagating cluster may either correspond to two or 
zero propagating cluster boundaries. The latter case will occur when the cluster is allowed 
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We refer to Refs. [29, 30] for more details. 
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to have cross-topology [43], i.e., to percolate across both periodic lattice directions. The 
corresponding bulk operator is denoted 0\ and has conformal weight 



hi(p) 



-3p 2 + 8p-A 
16p 



(5.13) 



Indeed, by duality this is nothing else than the order parameter (magnetization) operator. 

In the probabilistic regime (i.e., v > 0) the identity operator X is the most relevant (all 
other conformal weights are strictly positive). Similarly, one has h\ < hi, translating the 
fact that a single cluster is indeed most likely to have cross-topology. In the BK phase, 
v < and some of these statements need no longer be true, due to the presence of negative 
conformal weights. In fact, for any p > 2 the weight hi is the most negative, and so one 
would naively expect that the most "likely" configurations are those with two propagating 
cluster boundaries. This is however not true, due to a subtle effect: the amplitudes of the 
corresponding transfer matrix eigenvalues vanish identically on the torus [30]! 

So the remaining candidates for the most relevant operators are 0\ and Oi- From 
Eqs. (5.12)-(5.13) it is easy to see that the former (resp. the latter) will dominate for 2 < 
p < 4 (resp. for 4 < p < oo). This argument is of course valid in the thermodynamic 
limit only, but the conclusion will in fact hold true in finite size due to the underlying 
quantum group symmetry of the model. Accordingly we indeed observe numerically, for any 
width L and for both lattices, that branches of the limiting curves B cross the real g-axis 
at q = (p = 2) and q = 2 {p = 4); the former (resp. the latter) crossing corresponds to 
the equimodularity of the dominant eigenvalues of the t = and £ = 1 (resp. the £ = 1 and 
£ = 2) sectors. However, the phase containing q = 2 + does not extend to q = 4 (p = oo) 
for either lattice, since the chromatic line leaves the BK phase at go — 3 (square lattice), 
respectively q<i ~ 3.8196717312 (triangular lattice). 

We can further check the validity of this scenario by turning to the amplitudes of the 
dominant eigenvalues. For the order parameter operator 0\ the amplitude is q — 1 (by a 
simple counting argument), and this is indeed the amplitude of the dominant eigenvalue 
observed numerically for p e (2,4) for both lattices. (See Conjectures 5.3(b) and 5.8(b).) 
For the two-cluster operator (9 2 the amplitude should be either bf^ or [cf. Eqs. (2.32a)- 
(2.32b)]. And indeed the amplitude of the eigenvalues which is dominant in the phase 
containing q = 2 + is observed numerically to be b\ . (See Conjectures 5.3(c) and 5.8(c).) 
As a last check, note that from Eqs. (5.12)— (5.13) the effective central charge should read 



We have checked numerically that this is true on both lattices (using strips of width L < 12) 
to a very good precision; this check is essentially the contents of Figure 12(b). 

A final elusive point is the value of c in the interval q G (q 2 , qo) on the triangular lattice. 
The go(tri) = 4 case is known to be equivalent to Baxter's three-colouring problem on the 
hexagonal lattice and hence have c = 2. Our numerical data seem to indicate that the rest 
of the interval (q 2 , qo) may also correspond to a critical theory, the properties of which are 
unrelated to those of the BK phase. 




for 2 < p < 4 



for 4 < p < p 



(5.14) 
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C (l) = 2 - ~2 no ( 5 - 16 ) 



We have made an attempt to obtain an expression for the central charge in the interval 
q G (q2, qo) by performing some fits to our free-energy data. Our first two-parameter Ansatz 
was given by Eq. (5.6) without any correction-to-scaling term. In Figure 13 we show our 
estimates as a function of q for different values of L min = 6, 9, 12. Much better estimates 
can be obtained by including higher-order corrections to the above Ansatz. In previous 
work [9, 10, 13], we have found that including a correction term ~ L~ A greatly improves the 
estimates for the central charge. Thus, we have used the improved three-parameter Ansatz 

h(q) = /bu. k (?) + ^ + Ar 4 (5.15) 

In Figure 13 we have also displayed the estimates for c(q) obtained with this Ansatz and 
£ m in = 6, 9. Their stability is superior to that of the two-parameter estimates. The solid 
line represents a guess for the exact value of the (effective) central charge in this regime (see 
(5.19) below) 

28 

p 2 -~2p~- 12 

where p G (2, oo) is related to q by (5.10). 

Let us now explain in detail how we have obtained the expression (5.16). The functional 
form of the Ansatz (5.16) is motivated by already known results of the central charge in other 
regimes [c.f., (5.11)]. Thus, our goal is to obtain the central charge c(q) as a rational function 
of p with integer coefficients. We started with the estimates coming from the three-parameter 
Ansatz (5.15). It is more useful to plot the results as a function of 1/p, so that q = 4 (p = oo) 
corresponds to the origin. In Figure 14(a) we display the difference c(q = 4) — c(q) versus 
1/p for the two available data sets L min = 6, 9. Both data sets fall approximately on a single 
curve, especially on the interval corresponding to q G (q2,Qo)- This empirical result implies 
that the largest FSS correction to our estimates for c(q) corresponds to a global shift. We 
then fit the data for c(q = 4) — c(q) to a polynomial Ansatz. Our best result is given by 

c(q = 4)-c(q) = 28.60506^ + 32.11981^- + 897.88358^- (5.17) 

p2 pi p4 

This Ansatz is motivated by several facts: by definition, the curve should go through the 
origin, and it is clear in Figure 14(a) that its slope at 1/p = is zero. Furthermore, the 
addition of more terms to (5.17) only leads to less stable estimates. As our goal is to obtain 
a rational function of p, we computed the Pade approximant [2, 2] of the above polynomial 
(5.17). The result is given by 

t a\ i \ 28.60506 
° iq = A) - ° iq) = p2- 1.12287p- 30.12813 (5 ' 18) 

Finally, as we expect that the coefficients of the above rational function are integers, we 
systematically searched the set of closest integers which best approximates the data of Fig- 
ure 14(a). Our best candidates were 

c(? = 4)-c(9) = p 2_f p _ 12 (5- 19a) 

29 

c{q = A)-c{q) = p2 _ p _ 20 (5.19b) 
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In Figure 14(a) it is clear that both expressions give an excellent fit to the whole set of data 
points. In order to tell which integer-coefficient rational approximant was best, we plotted 
the difference between our estimates and the fits. In Figure 14(b) we show this difference for 
the above approximants (5.19). We observe that the approximant (5.19a) is slightly superior 
to the other one (5.19b). 

The final step is to fix the overall shift, i.e., the value of c(q = 4). As the triangular-lattice 
Potts antiferromagnet at zero temperature is critical and can be mapped onto a Gaussian 
model with two bosons [44,45], we expect that 



Thus, we arrive at the final expression (5.16). Indeed, this computation should be taken with 
a grain of salt: we have found a function with the expected properties (rational function of 
p with integer coefficients) that describes rather accurately our estimates for the central 
charge of our model on the interval q G (q2,Qo)- This does not mean that this is unique 
possible solution with those properties nor that the right expression for c(q) should satisfy 
those properties. 15 Another important point is that our conjectured expression for the central 
charge, although it was first devised to deal with the interval q G (g 2 , qo), is able to reproduce 
rather accurately the data outside the latter interval. 



At the end of Section 5.3, we have already discussed the branch of the triangular-lattice 
limiting curve El going through the point q — 2, and separating the i = 1 and I = 2 phases. 
In this subsection we will study further properties of this branch in both the square and 
triangular limiting curves. 

The first question we want to address is that of the density of chromatic zeros along this 
branch for finite strips of size Lx N. In order to estimate this density, we first fixed a region 
in the complex g-plane containing the q = 2 branch. This region is 



In Figure 15, we show the density of the chromatic zeros p in the regions (5.21) for some 
finite strips of size L x iV as a function of the aspect ratio L/N. At fixed aspect ratio, we 
expect the density p to converge to some (small) positive constant. In Figure 15, we observe 
that for fixed L/N, the density is decreasing in L, and it attains at L = 7 very small values 
< 0.02. If this monotonicity holds for larger values of L, we expect that the density of zeros 
along the q = 2 branch will converge to zero in the thermodynamic limit (in contrast with 
the behavior of the density of zeros in other parts of the limiting curve). 

Another way to tackle this problem is consider the parameter t introduced in Ref. [11, 
Section 4.1.1]. On the limiting curve Bl, there are generically two equimodular leading 
eigenvalues Ai and A 2 . Indeed, their ratio is just a phase 



c(q = 4) = 2 



(5.20) 



5.5 The "q = 2 branch" 



Sq =>. 1.8 < Reg < 2.2, 
Tri 1.8 < Reg < 2.6, 



Img| < 1.2 
Img| < 1.4 



(5.21a) 
(5.21b) 




8 



6eR 



(5.22) 



15 Accordingly, we have referred to (5.16) as a guess rather than a conjecture. 
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In Ref. [11], instead of 9, we found it more convenient to work with the real parameter t, 
defined as 

t = tan ^0 , t e [0,oo) (5.23) 

In practice, we use this parameter to ensure the plots of the limiting curves Bl are correct 
and we do not miss any small feature. Note that t = 9 = corresponds to two identical 
eigenvalues. At q — 2, we have that t — for any L > 3; but as we move along the 
q = 2 branch, t increases. Let us now study how the parameter t evolves (along this 
branch) as a function of \q — 2|. In Figure 16, we show (for both the square and triangular 
lattices) the values of t along the q = 2 branch for points in the complex g-plane with 
Img = 0, 0.01, . . . , 0.10. It is clear that for fixed L, the relation is approximately linear. For 
the square lattice, the slope is given approximately by 1/L, while for the triangular lattice 
the slope is given approximately by \/3/(2L) + 0.653/L 3 . In both cases, the discrepancies 
between the estimated slopes and the above results are of order O(10 -4 ). Thus, we expect 
that in the thermodynamic limit both leading eigenvalues Ai and A2 become equal along this 
branch. 

Finally, we can study the behavior at q = 2 by considering the series expansion of the two 
equimodular leading eigenvalues Ai and A2 (corresponding to the £ = 1 and £ = 2 sectors, 
respectively). Indeed, such series expansions are expected to exist, as the eigenvalues \(q) 
are algebraic functions of q. On each phase, we can compute the series expansion of the 
leading eigenvalue at q = 2 + e 

00 

A,(g + e) = $>,, fc e fc (5.24) 

k=0 

For L = 2,3, this can be done symbolically by using mathematica. For 4 < L < 7, we can 
only obtain (again with mathematica) the series numerically, as the symbolic expressions for 
the coefficients are too lengthy. In practice, we are able to obtain the first three coefficients 
Aj,fc (k = 0, 1,2) in (5. 24). 16 With these coefficients, we can obtain the corresponding series 
expansion for the free energy: 

U log 1^(2 + 6)1 
h(2 + e) = \ 

|^log|A 2 (2 + e)| 

For 4 < L < 12, we numerically obtained the contribution to the free energy /l,i(?) = 
(1/L) log |Aj(g)| of the two leading eigenvalues Aj close to q — 2. In particular, we obtained 
the values fi,i{q = 2) and fi,i{q = 2± 10 -5 ) for % — 1, 2. From these six values, we estimated 
the first three coefficients (with k = 0,1,2) in (5.25). Indeed, we know that the free 

energy is continuous at q — 2, thus we expect that A/l(2) = /| 2 q — f^ = 0, as we find in the 
actual computation. We are interested in computing the increment of the first derivatives 
of the free energy /l(?) at q = 2. In practice, we have computed the increments of the first 



For the triangular lattice, this computation could only be performed up to width L = 6. 



E/Se fc if e < 

k=0 

tf ( L^ k if e > 

fc=0 



(5.25) 
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two derivatives: 17 

El 
Cj 




q=2 y 



d'fL 



dq 1 



/S " /$ (5.26a) 
= 2(iS - /S) (5.26b) 



5=2, 



If £x 7^ 0, then the transition at q = 2 is of first order. However, if E L = and C L ^ 0, then 
the transition is of second order (and this argument can be generalized to higher-order phase 
transitions). In Figure 17 we have plotted the coefficients El and Cl versus 1/L for both 
the square and triangular lattices. The plots suggest that both coefficients for both lattices 
tend to a limit as L — > oo. If we fit the S^-data to a power law El = -Eoo + AeL~ Ue , we 
find that l-E^I ~ 10~ 4 is indeed very small (in modulus) for L min = 10: 

Sq £ L « 9.8 x 10~ 7 - 1.98329 L' 1 " 737 (5.27a) 
Tri E L « -0.00010 - 2.12401 L~ 2 - 02506 (5.27b) 

In particular, our square-lattice data can be fitted very well with the simple Ansatz: £x(sq) = 
—2L~ 2 . These fits are depicted in the corresponding plots in Figure 17. We have chosen 
L min to be the largest possible value, so that there are no effective degrees of freedom in the 
fits. 

The power-law fits Cl = C^ + AqL^ are not as stable as the preceding ones: we find 
sizeable parity effects in the estimates, so we have to split the whole data set into two smaller 
sets with well-defined L-parity. Thus, for each value of L min , we have fitted the data with 
L = L min , L min + 2, L min + 4. The estimates for even L are given for L min = 8 by: 

Sq C L w 0.00051 - 0.35387 L- 1 - 73364 (5.28a) 
Tri C L ~ -0.00012 - 0.99690 L" 2 18046 (5.28b) 

For odd L with L min = 7, the corresponding estimates are: 

Sq C L ~ 0.00116 - 0.30522 L~ L62941 (5.29a) 
Tri =>■ C L w -0.00092 - 1.42076 L" 2 - 38725 (5.29b) 

Our preferred fits are those given by (5.28), as L min is larger; these fits are depicted in the 
corresponding plots in Figure 17. They also show show that C L tend to some small (in 
modulus) number in the thermodynamic limit |Coo| ^5x 10 -4 . We conjecture that for both 
lattices, El,Cl — > as L — > oo. Indeed, this does not mean that there is no transition at 
all at q = 2 in the thermodynamic limit, but that the transition (if any) is at least of third 
order. 

In this section we have presented three numerical approaches to study the thermodynamic 
limit along the q = 2 branch. None of them is conclusive; but taken together they give strong 



17 For 4 < L < 7, we had two distinct ways of estimating the coefficients El and Cl, namely symbolically 
and numerically. We checked that the difference between those estimates was < 1CP 11 for El, and < 1CP 6 
for C L . 
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support to the following scenario: even though at finite L there are two different phases 
coexisting along the q = 2 branch, in the thermodynamic limit there is a unique phase. This 
scenario is consistent with the three results discussed above: in the thermodynamic limit 
L — > oo: 1) The density of chromatic zeros goes to zero along this branch; 2) The ratio of 
eigenvalues along this branch tend to one A2/A1 — > 1; and 3) At q = 2 the increment in the 
first two derivatives of the free energy fi{q) (with respect to q) tends to zero. Thus, this 
branch seems to be somehow an artifact of the boundary conditions; in particular, of the 
periodic longitudinal boundary conditions, as this branch also appears in the cyclic case [9], 
and it is absent for free and cylindrical boundary conditions [11-13]. 18 This conclusion agrees 
with Baxter's solution [14] in the sense that only three phases survive in the phase diagram 
of the triangular-lattice chromatic polynomial. 

6 Conclusions 

In this paper we have extended our transfer-matrix studies [9,11-13] to toroidal boundary 
conditions. We have obtained the exact expression of the chromatic polynomial for square- 
and triangular-lattice strips of fixed width 2 < L < 7 and arbitrary length N. The cases 
L = 5, 6, 7 are new in the literature (note that going to larger widths L > 8 would require 
better computer facilities and/or improved symbolic algorithms). We have also obtained the 
limiting curves Bl and isolated limiting points for these strips in the infinite-length limit 
N — > 00. The limiting curve Bl for a fixed value of L provides the exact phase diagram for 
this semi-infinite strip. 

By studying the features of these finite-L phase diagrams and using techniques from 
CFT and FSS, we have been able to determine some properties of the chromatic-polynomial 
phase diagram in the thermodynamic limit L — > 00. These findings are summarized in the 
following points: 

• Toroidal boundary conditions appear to give results qualitatively and quantitatively 
closer to those of the thermodynamic limit than the other standard boundary condi- 
tions (free, cylindrical, and cyclic) [9,11-13]. In particular, the exact value of q c is 
obtained for all square lattices of (finite) odd widths, and for all triangular lattices of 
(finite) widths not a multiple of three. 

• As for cyclic boundary conditions [9], we find that the limiting curves Bl always enclose 
regions of the complex g-plane, thus determining true phases as in the expected infinite- 
volume phase diagram. Each phase is also characterized by the number of bridges t 

18 A similar (but not related) phenomenon appears in a recent study of the RSOS [8] representation of 
the Potts model when q = B p with p = B p (1.8). In the complex u-plane, the large-|v| region (which 
corresponds to the high-temperature regime of the model and is expected to be non-critical) displays several 
non-bounded outward branches separating different phases. The origin of these branches is simple: there 
are two almost-degenerate eigenvalues, each of them corresponding to a different topological sector. These 
two eigenvalues become exactly equimodular along those non-bounded outward branches, whose number 
increases with L. While for any finite L there are several phase-transition lines in the large- 1> region, in the 
thermodynamic limit we expect that both eigenvalues become exactly identical over all this region. Thus, we 
recover the standard result that there are no phase-transition lines in the large-u (high-temperature) regime 
of this model. 
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associated to the leading eigenvalue; but for toroidal boundary conditions only the 
i = 0, 1,2 sectors contribute for both lattices. (For the triangular lattice with cyclic 
boundary conditions we expect to have contributions up to I = 6 bridges). The other 
difference with the cyclic triangular strip is that there are now two different phases 
characterized by £ = 2 bridges (For cyclic boundary conditions, there is at most one 
phase for each value of £). 

It is rather remarkable that for cyclic boundary conditions [9,29], all the Beraha num- 
bers Bp (1.8) played a special role in the BK phase, either as phase transition points 
(even p) or as isolated limiting points (odd p). In contrast, for toroidal boundary 
conditions only the finite subset of integer q play such special roles: there are phase 
transitions at q = 0, 2, 4 and isolated limiting points at q — 1, 3. This observation alone 
is a nice illustration of the crucial influence of boundary conditions on antiferromag- 
netic systems. Note however, that we still expect (and have to some extent numerically 
observed) that the full set of B p will be the loci of level crossings of sub-dominant eigen- 
values. 

Seeing that the maximum chromatic number H{g) of a graph G increases with the genus 
g of the surface that it is embedded in according to (1.5), one might have expected 
that the toroidal topology studied here would have revealed further structure of the 
limiting curves for q G (4,7). This is not the case. The reason is that it is precisely 
the interplay between the properties of the BK phase and the representation theory of 
the quantum group U t (SU(2)) for t a root of unity that is responsible for generating 
the partition-function zeros. This mechanism is confined to the interval q G [0,4]. The 
graphs having chromatic zeros for q G (4, 7) must therefore be rather exceptional and 
have little to do with the recursive families of graphs studied here. One might even 
speculate that the number of such exceptional graphs might be small in some sense, 
or even finite. These conclusions are expected to carry over to higher genus as well. 

For toroidal triangular strips we have found two phase transitions inside the critical 
region: q = 2 and q = q^ (1.11). 

This result is rather puzzling when compared to Baxter's exact solution [14] and the 
re-analysis discussed in Ref. [9] of Baxter's eigenvalues. On one hand, the value of ge 
coincides with that predicted by Baxter based on the result that only three eigenvalues 
are relevant in the complex g-plane. This conclusion stems from the Bethe-Ansatz 
solution of a closely related model [14]. The mapping between these two models is 
exact, except for the boundary conditions. Usually boundary conditions have no effect 
in the ferromagnetic regime; but here we are deep in the antiferromagnetic one. 

On the other hand, a careful and high-precision analysis of these three eigenvalues [9] 
revealed that, although the outer boundary of Baxter's infinite-volume curve £>oo was 
indeed correct, the inner structure was richer: there were infinitely many more branches 
than those included by Baxter and those indicated by the finite-strip chromatic zeros. 
In particular, the point ge (1-H) does not belong to the limiting curve £>oo, and its role 
is now played by qc (1.12). 

In Section 5.5 we have found some numerical evidence that the branch of the limiting 
curves Bl going through q = 2 might disappear in the true thermodynamic limit 
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L — > oo. This may provide an explanation of why this branch is absent from the 
triangular-lattice limiting curve B^ proposed by Baxter [14]. However, it is still not 
clear how to explain the existence of the q = 2 branch for finite widths from Baxter's 
analysis. The resolution of these discrepancies would most likely involve identifying 
(at least) one more analytic expression for the free energy from the Bethe Ansatz 
equations of Ref. [14]. Such a study might also shed some light on the possible existence 
of additional enclosed regions in B^, such as the triangular-shaped regions discussed 
towards the end of Section 5.2. 

• In the triangular-lattice case we obtain a novel behavior: one important feature of 
the infinite- volume phase diagram does depend on the parity of the width (i.e., on 
L mod 3). Our numerics strongly suggests that 



In the square-lattice case, we have that this limit for even widths is equal to <? c (sq) = 
q (L = 2p+l) for all p > 1. 

• In the light of the numerical results given at the end of Section 5, we might finally 
speculate that the full solution of the triangular-lattice chromatic polynomial (in the 
sense of a Bethe Ansatz solution) should reveal not only two distinct bulk free energies 
for q E [0,4], but also two different sets of associated critical scaling levels throughout 
that interval. In other words, we suggest that there exists two "superposed" critical 
theories throughout the interval, one having the physics of the BK phase, and the 
other being some deformation of the two-boson theory [44, 45] responsible for c = 2 
at q — 4. This is a strong motivation for a more refined study of the Bethe equations 
corresponding to Baxter's g 2 solution [14]. Note that these expectations apply also to 
the Nienhuis loop model [41]. 

A Two combinatorial identities 

In this section we prove the following combinatorial identities that were used in Sec- 
tion 2.2. 

Lemma A.l For all integer m > 1, the following identities hold 




if L mod 3^0 



if L mod 3 = 



(6.1) 




(A.lb) 



(A.la) 
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Proof. 

Let us start with the proof of (A. la). The first goal is to make explicit the m factor in 
the r.h.s. of (A. la) 



J2p 



2m 
m — p 



(2m 



p=0 



m 



E 

,p=0 



2m 
P 



m 

p=0 



2m 
P 



- 2 E 

p=0 

2m- 1 
P 



2m\ ( 2m — 1 



P 



p — 1 



2m - 1 
p-l 



(A.2) 



where we have first changed the summation index p —> m — p, then used the absorption 
identity [46, Eq. (5.5)] 





a 








= i 
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and finally used the addition formula [46, Eq. (5.8)] 
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If we use (A.4) to expand ( 2m ) in (A.2), we arrive at the final result: 
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Let us now prove Eq. (A. lb). The idea is similar to the above proof; but now it is more 
efficient to extract a common factor (m + 1) rather than m. Thus, if we follow the same 
steps as above, we get 
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We now use (A. 4) to expand the two binomials inside the brackets in (A. 6), and we get the 
desired result: 




where the last equality follows from the fact that m is a positive integer. | 
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Table 1: Dimensionality of the transfer matrix for toroidal boundary conditions. For each 
width m, we give the number of non-crossing connectivities C m)Tn (2.17), the number of 
non-crossing connectivities C m+ \, m (2.19), the number of non-crossing no n- nearest-neighbor 
connectivities -D m m , the number of generated connectivity classes modulo reflections (resp. 
reflections and translations) of non-crossing non-nearest-neighbor connectivities TriTorus' 
(resp. SqTorus'), and the number of distinct eigenvalues for a triangular-lattice (resp. square- 
lattice) strip with toroidal boundary conditions TriTorus (resp. SqTorus). 
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Table 2: Block structure of the transfer matrix of a toroidal square-lattice strip of width m 
as a function of number of bridges t. For each strip width m, we quote the dimension of 
the diagonal block T^(mx) for a given number of bridges £, and the total dimension of the 
transfer matrix SqTorus'(m) = dimT(m T ). 
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Table 3: Block structure of the transfer matrix of a toroidal triangular-lattice strip of width 
m as a function of number of bridges i. For each strip width m, we quote the dimension Ng 
of the diagonal block Ty(mT) for a given number of bridges £, and the total dimension of 
the transfer matrix SqTorus'(m) = dimT(m T ). 
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Table 4: Number of distinct eigenvalues for a toroidal square-lattice strip of width m as a 
function of number of bridges £. For each strip width m, we quote the number of different 
eigenvalues L$ for a given number of bridges £, and the total number of distinct eigenvalues 
SqTorus(m). Notice that the sum YlT=o^^ m ) ma y ^ e g rea + er than SqTorus(m), as a few 
eigenvalues can appear on different £ blocks. The results marked with a t are conjectures 
based on the numerical evaluation of the corresponding eigenvalues. 
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Table 5: Number of distinct eigenvalues for a toroidal triangular-lattice strip of width m as 
a function of number of bridges £. For each strip width m, we quote the number of different 
eigenvalues Le for a given number of bridges £, and the total number of distinct eigenvalues 
TriTorus(m). Notice that the sum X^o-^( m ) ma y ^ e g rea + er than TriTorus(m), as a few 
eigenvalues can appear on different £ blocks. The results marked with a ' are conjectures 
based on the numerical evaluation of the corresponding eigenvalues. 
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Table 6: Summary of qualitative results for the square-lattice eigenvalue-crossing curves B 
and for the isolated limiting points of zeros. For each toroidal square-lattice strip considered 
in this paper, we give the number of connected components of B C), the number of 
endpoints E), the number of T points (# T), the number of double points D), and 
the number of enclosed regions (# ER). We also give the value qo(L) which is the largest 
real value where B intersects the real axis, and the number of real isolated limiting points of 
zeros (# RI). 
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Table 7: Summary of qualitative results for the triangular-lattice eigenvalue-crossing curves B 
and for the isolated limiting points of zeros. For each toroidal square-lattice strip considered 
in this paper, we give the number of connected components of B (# C), the number of 
endpoints (# E), the number of T points (# T), the number of double points (# D), and 
the number of enclosed regions ER). We also give the value qo(L) which is the largest 
real value where B intersects the real axis, the value q2(L) where there is a crossing between 
two dominant eigenvalues that belong to the £ = 2 sector, and the number of real isolated 
limiting points of zeros (# RI). 
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Table 8: Summary of fits of the go(3p) and g2(3p) data (1 < p < 4) for the zero-temperature 
triangular-lattice Potts antiferromagnet. For each data set and Ansatz, we show the different 
estimates as a function of the parameter L min . See text for an explanation of how the fits 
were performed. 
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Figure 1: The complete graph K 7 embedded on a torus. The solid lines represent the edges 
of the graph. To build the torus, we should glue together opposite dashed lines. We thank 
David Wilson for helping us preparing this plot. 
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Figure 2: Square lattice with toroidal boundary conditions of size 4x2. This lattice is 
obtained from a square lattice with cyclic boundary conditions of size 4 x 3 by identifying 
the top and bottom rows i <-> i' (i — 1, . . . , 4 = m). 




Figure 3: Triangular lattice with toroidal boundary conditions of size 4x2. This lattice 
is obtained from a triangular lattice with free boundary conditions of size 5 x 3 by first 
identifying the first and last columns (i.e., 1 5 and V <-> 5'), so we get a triangular lattice 
of size 4x3 with cylindrical boundary conditions. Then, as for the square lattice, we identify 
the top and bottom rows i <-> i' (i — 1, . . . , 4 = m). 
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(b) 



Figure 4: Connectivity state for a toroidal strip of width 4. (a) We show the connectivity 
state V = 3' 52,4,4', which can be seen as a bottom-row connectivity state V^ot — ^i',3' 
(lower solid thick line), a top-row connectivity state V t0 p = #2,4 (upper solid thick line), and 
a bridge connecting both rows 64^ (dashed thick line). The thin solid lines show the top 
(resp. bottom) row containing the unprimed (resp. primed) sites. The thick arrow shows the 
transfer direction (upwards), (b) The same connectivity state depicted in such a way that 
the transverse periodic boundary conditions are clearer. The transfer direction is outwards 
(thick arrow). 
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Zeros sq lattice L = 2 T Zeros sq lattice L = 3 T 




-10 12 3 -10 12 3 

Re(q) Re(q) 

(c) (d) 

Figure 5: Limiting curves for square-lattice strips of width (a) L — 2, (b) L — 3, (c) L = 4, 
and (d) L — 5 with toroidal boundary conditions. We also show the zeros for the strips 
L-p x (5L)p (black □) and L-p x (10L)p (red o) for the same values of L. 
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Zeros sq lattice L = 6 T 



Zeros sq lattice L = 7 T 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

-10 12 3 -10 12 3 

Re(q) Re(q) 

(a) (b) 

Figure 6: Limiting curves for square-lattice strips of width (a) L — 6, and (b) L — 7 with 
toroidal boundary conditions. We also show the zeros for the strips Lp x (5L)p (black □) 
and L F x (10L) P (red o) for the same values of L. 
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Zeros tri lattice L = 2 T Zeros tri lattice L = 3 T 




Re(q) Re(q) 

(c) (d) 

Figure 7: Limiting curves for triangular-lattice strips of width (a) L — 2, (b) L — 3, (c) 
L — 4, and (d) L — 5 with toroidal boundary conditions. We also show the zeros for the 
strips L-p x (5L)p (black □) and Lp x (10L)p (red o) for the same values of L. 
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Figure 8: Limiting curves for triangular-lattice strips of width (a) L = 6, and (b) L = 7 with 
toroidal boundary conditions. We also show the zeros for the strips Lp x (5L)p (black □) 
and Lp x (10L)p (red o) for the same values of L. 
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Limiting Curves Sq Lattice 



t — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — r 




j i i i I i i i i I i i i i I i i i i I L 



-10 1 2 3 

Re(q) 

Figure 9: Limiting curves for square-lattice strips of width L and toroidal boundary condi- 
tions. We show the results for L = 2 (black), L = 3 (red), L = 4 (green), L = 5 (blue), 
L = 6 (light blue), and L = 7 (pink). We also show the limiting curve for a square-lattice 
strip of width L = 10 with cylindrical boundary conditions (orange) [12]. 
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Limiting Curves Tri Lattice 




Figure 10: Limiting curves for triangular-lattice strips of width L and toroidal boundary 
conditions. We show the results for L = 2 (black), L = 3 (red), L = 4 (green), L = 5 (blue), 
L = 6 (light blue), and L = 7 (pink). We also show the infinite-width limiting curve 
obtained by Baxter [14] (orange). 
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q 

Figure 11: Free energy per unit area fh{q) for the triangular-lattice chromatic polynomial 
on an infinitely long strip of width L — 12 and toroidal boundary conditions. We show the 
free energy associated to a selected set of eigenvalues of the transfer matrix T(12 T ) as 

a function of q G (3.5,4). For each sector £ = 0, 1,2, we show the dominant and the first 
sub-dominant eigenvalues respectively with colors red/orange for i = 0, blue/cyan for £ = 1, 
and black/dark gray for i = 3. Those eigenvalues that are decreasing (resp. increasing) in q 
are depicted with a solid (resp. dot-dashed) curve. For comparison, we also show Baxter's 
free energies (1/G) log #2 and (1/G) log gs [14] normalized per unit area and depicted as solid 
pink and violet, respectively. The vertical dashed dark gray line shows Baxter's value ge 
(1.11). 
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Figure 12: Ratios of the free energy per unit area for an infinitely long triangular-lattice 
strip of width L and toroidal boundary conditions and Baxter's free energy per unit area 
(1/G) log 03 (g) [14]. (a) Ratio of the free energy per unit area /z,(?) (5.4). (b) Ratio of the 
corrected free energy per unit area / COT r,i(9) (5-7). For each strip width L = 6 (orange), L = 9 
(pink), L = 12 (blue), and L — 15 (black), we show the ratio of the free energies associated 
to the first three levels £ = (solid line), £ = 1 (dashed line), and £ = 2 (dot-dashed lines). 
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Central charge tri lattice 




Figure 13: Central charge for the triangular-lattice chromatic polynomial on infinitely long 
strips with toroidal boundary conditions. We show several two-parameter fits (depicted in 
red): L = 6, 9 (+), L = 9, 12 (*), and L = 12, 15 (x). We also show three-parameter fits 
using the Ansatz (5.15) (depicted in black): L = 6,9,12 (o), and L = 9,12,15 (□). The 
solid blue line represents our guess for the central charge (5.16). The vertical dashed dark 
gray line shows Baxter's value ge (1-11). 
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Figure 14: Central charge for the triangular-lattice chromatic polynomial on infinitely long 
strips with toroidal boundary conditions, (a) Difference c(q = 4) — c(q) as a function of 
1/p [c.f., (5.10)]. We show two estimates coming from the three-parameter fits (5.15) with 
L = 6,9, 12 (black o) and L = 9, 12, 15 (red □). The (almost coincident) solid lines repre- 
sent our two best candidates for the central charge (5.19a) (blue), and (5.19b) (red), (b) 
Error between our numerical estimates for the conformal charge and our conjectured expres- 
sions (5.16). We show the difference c(q) — cgt(g) between our three-parameter fits and our 
conjectured expressions Cfi t (g). Points in black (resp. red) are obtained from the fits with 
L = 6,9,12 (resp. L = 9,12,15). Solid [resp. empty] symbols are obtained with Cfi t given 
by (5.19a) [resp. (5.19b)]. In both plots, the vertical dashed dark gray line shows Baxter's 
value 1/pb w 0.0681036688 obtained from (1.11)/(5.10). 
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tri lattice: q=2 branch 
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Figure 15: Density of chromatic zeros for (a) square- and (b) triangular-lattice strips with 
toroidal boundary conditions along the q = 2 branch of the limiting curve B^. In order to 
compute the density of zeros, We have considered the regions given by (5.21). For a given 
lattice strip of width L and length N, we plot the ratio of zeros falling in the respective 
region as a function of the inverse aspect ratio (i.e., L/N). The widths are L = 2 (black □), 
L = 3 (red o), L = 4 (violet A), L = 5 (blue o), L = 6 (orange A), and L = 7 (pink ♦). 
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Figure 16: Values of the parameter t (5.23) along the q = 2 branch of the limiting curve Bl 
for the (a) square- and (b) triangular-lattice strip. For a lattice strip of width L, we show 
the value of t as a function of \q — 2|. The widths are L = 3 (red o), L = 4 (violet A), L = 5 
(blue o), L = 6 (orange A), and L = 7 (pink The slopes of the depicted lines are given 
approximately by 1/L (square) and y/3/(2L) +0.653/L 3 (triangular). 
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Figure 17: Coefficients of the series expansions for the increment of the first (a) and second 
(b) derivatives of the free energy with respect to q evaluated at q — 2. In each plot we 
show the corresponding coefficient for the square (red o) and triangular (blue A) lattices as 
a function of 1/L. The solid lines depict the best power-law fits to the corresponding data 
set (5.27)/(5.28). 
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